arXiv: 1506.04662vl [math.OC] 15 Jun 2015 


OPTIMAL CONTROL OF THE SWEEPING PROCESS 
OVER POLYHEDRAL CONTROLLED SETS 
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Abstract. The paper addresses a new class of optimal control problems governed by the dissipative and dis¬ 
continuous differential inclusion of the sweeping/Moreau process while using controls to determine the best shape 
of moving convex polyhedra in order to optimize the given Bolza-type functional, which depends on control and 
state variables as well as their velocities. Besides the highly non-Lipschitzian nature of the unbounded differential 
inclusion of the controlled sweeping process, the optimal control problems under consideration contain intrinsic 
state constraints of the inequality and equality types. All of this creates serious challenges for deriving necessary 
optimality conditions. We develop here the method of discrete approximations and combine it with advanced tools 
of first-order and second-order variational analysis and generalized differentiation. This approach allows us to estab¬ 
lish constructive necessary optimality conditions for local minimizers of the controlled sweeping process expressed 
entirely in terms of the problem data under fairly unrestrictive assumptions. As a by-product of the developed 
approach, we prove the strong VT^’^-convergence of optimal solutions of discrete approximations to a given local 
minimizer of the continuous-time system and derive necessary optimality conditions for the discrete counterparts. 
The established necessary optimality conditions for the sweeping process are illustrated by several examples. 
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1 Introduction and Problem Formulation 

This paper is devoted to the study of the sweeping process, a class of models introduced by Jean-Jacques 
Moreau in the 1970s to describe a number of quasistatic mechanical problems; see [28, 30, 31] and the book 
[22] for more details. Besides the original motivations, models of this type have found significant applica¬ 
tions to elastoplasticity [15], hysteresis [18], electric circuits [1], etc. For its own sake, the sweeping process 
theory has become an important area of nonlinear and variational analysis with numerous mathematical 
achievements and challenging open questions; see, e.g., [13, 20] and the references therein. 

Mathematically the sweeping process is governed by the dissipative differential inclusion 

(1.1) x(t) € —N[x{t);C{t)) a.e. fG[0,T], 

which describes the movement of a point belonging to a continuous moving set C(t) while its velocity 
belongs for a.e. t to the negative normal cone to C{t) at x{t). The Cauchy problem a;(0) = xq for the 
sweeping process (1.1) enjoys a developed well-posedness theory for convex and mildly nonconvex moving 
sets; see, e.g., [13]. Higher-order and state-dependent sweeping processes have also been studied (to a 
much lesser extent) in the literature; see, e.g., the book [9]. Let us also mention the recent paper [21], 
which contains existence and well-posedness results obtained via advanced tools of variational analysis and 
generalized differentiation for a broad class of evolution systems including the sweeping process. 
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Among the central issues of the sweeping process theory is establishing the existence and uniqueness of 
solutions to the Cauchy problem for the sweeping differential inclusion (1.1) under reasonable assumptions 
on the given moving set C{t). This tells us that it does not make any sense to optimize the sweeping process 
generated by the given set C(t) from the standard viewpoint of optimal control theory well developed for 
Lipschitzian differential inclusions and the like; see, e.g., [25, 38, 42] and the references therein. 

In our first paper on the sweeping process [11] we suggested to take a new viewpoint on optimizing the 
dynamical system (1.1) by controlling the moving set C(t) with the usage of control actions that change the 
shape of C{t) and hence the right-hand side of the sweeping differential inclusion (1.1). This idea was partly 
implemented in [11] for the case when the sweeping process was driven by a moving affine hyperplane whose 
normal direction and boundary were acting as control variables. Furthermore, it was assumed in [11] the 
independence of the running cost on time, control variables, and control velocities as well as the uniform 
Lipschitzian continuity of feasible controls. Apparently this first attempt was limited from both viewpoints 
of control theory and possible applications. A much more realistic while significantly more challenging case 
appears when, along with general running costs depending on state, control, and their velocity variables, 
controlled moving sets are described by convex polyhedra governed by finitely many controls in normal 
directions and polyhedron boundaries under inequality and equality constraints. Polyhedral descriptions 
of moving sets in the (uncontrolled) sweeping process were largely explored, e.g., in [15, 18, 19], where the 
reader can find interesting applications to particular models of elastoplasticity and hysteresis. 

In the other line of development we mention the recent paper [8] and its subsequent extension [2], which 
address a different class of optimal control problems for an equivalent variational inequality description of 
the sweeping process of the rate-independent hysteresis type, where the convex moving set is fixed while 
controls appear in an associated ordinary differential equation. Another recent paper [10], in the framework 
of BV solutions of a sweeping process whose given moving set is lower semicontinuous with nonempty 
interior, concerns relaxation issues and dynamic programming. Controls appear there via perturbations of 
the dynamics given by the normal cone while being the barycenter of a Borel finite measure. 

In this paper we study the following optimal control problem (P) of the generalized Bolza type for the 
sweeping process (1.1) as well as some of its modifications. Given an extended-real-valued terminal cost 
function K” —K := (— 00 , 00 ] and a running cost i: [0,T] x minimize the functional 

(1.2) J[x,u,b]-. = (p(^x{T)^ + f i(t,x{t),u{t),b{t),x{t),u{t),b{t)) dt 

Jo 

over the controlled sweeping dynamics described by 

(1.3) x{t) € —N[x{t);C{t)) for a.e. t G [0,T], x{0) := xq G C{0) 

with the inequality and equality constraint defined by 

(1.4) C{t) := {a; £ {ui{t), x) < bi(t), i = 1,..., m} 

(1.5) with j|Mi(t)jj = 1 for all t G [0,T], i = 1,... ,m, 

where the controls actions u{-) = (ui{-),... and b{-) = (bi{-),... ,bm{-)) are absolutely continuous 

on [0, T], the final time T is fixed, and the absolutely continuous trajectories x{-) of the differential inclusion 
are understood in the standard sense of Caratheodory. This class of problems contains several novel 
features, which either have never been investigated or have been studied insufficiently in control theory; 
see more discussions below. Recall now that the normal cone to a convex set C G1 H." is defined by 

(1.6) N{x-,C):={vGM.'^\{v,y-x)<0, y G C} ii x G C a,nd N{x-,C) := dl ii x ^ C. 

Hence the sweeping process with a given moving set (1.1) can be considered as an evolution variational 
inequality, or a differential variational inequality in the terminology of [32]. On the other hand, the control 
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model (1.3)-(1.5) relates rather to evolutionary quasi-variational inequalities with controlled parameters, 
which should be determined to optimize the dynamical process. 

The main goal of this paper is to derive necessary optimality conditions entirely in terms of the problem 
data for the so-called intermediate local minimizers of (P) and its modifications that occupy an intermediate 
position between the standard notions of weak and strong local minima in variational and control problems; 
see Section 3 for more discussions. Our approach is based on developing an appropriate version of the 
method of discrete approximations, which largely follows the scheme of [23, 25] implemented therein for the 
case of Lipschitzian and uniformly bounded differential inclusions, while now requiring a novel extension to 
the case of totally non-Lipschitzian and unbounded differential inclusions in (1.3). Some results on discrete 
approximations of feasible trajectories of (1.2) and the convergence of optimal solutions to appropriate 
discretizations of the continuous-time system in (1.2)-(1.5) have been recently obtained in our preceding 
paper [12]. However, they do not provide enough information for passing to the limit in necessary optimality 
conditions for discrete approximations and thus establishing in this way necessary optimality conditions 
for local minimizers of the original continuous-time systems governed by the controlled sweeping process. 

In this paper we are going to proceed further in this direction by improving the previous discrete 
approximation results to make it possible deriving necessary optimality conditions for the continuous-time 
systems by passing to the limit from those for their discrete approximations that are proved to satisfy the 
desired well-posedness and convergence properties. The realization of this approach requires overcoming 
significant difficulties, which have never been addressed earlier in control theory from this or any other 
method of deriving necessary optimality conditions even for more simple problems with smooth data. 
Besides the aforementioned totally non-Lipschitzian and unbounded nature of the sweeping process, serious 
challenges come, in particular, from the intrinsic presence of state constraints of the inequality and equality 
types combined with the quasi-variational inequality structure of the controlled sweeping process. Indeed, 
we show in Section 3 that problem (P) can be rewritten in the more conventional form of the generalized 
Bolza problem for a non-Lipschitzian and unbounded differential inclusion with a fixed right-hand side, 
where the relations in (1.4) and (1.5) are treated as state constraints of the inequality and equality type, 
respectively. It is worth mentioning that, in contrast to the inequality state constraints well studied for 
standard control systems and Lipschitzian differential inclusions (see, e.g., [4, 42] and the references therein), 
the equality state constraints have been just very recently addressed in [5] for smooth control systems under 
regularity assumptions formulated via full ranks of the corresponding Jacobians of the constraint functions. 
Needless to say that neither the results nor the approach of [5] can be applied in our setting. 

A crucial ingredient of our approach within the method of discrete approximations is applying advanced 
nonconvex tools of first-order and second-order of variational analysis and generalized differentiation, which 
are required even in the case of smooth terminal and running costs in the convex cone setting of (1.1). 
This allows us not only to establish the desired strong convergence of discrete approximations and then to 
pass to the liming in the necessary optimality conditions obtained for discrete problems, but also to derive 
necessary optimality conditions for the continuous-time control system entirely in terms of the problem 
data and the given local optimal solution to the controlled sweeping process. 

In fact, our major necessary optimality conditions are derived for a certain parametric perturbation 
{P'^) of the original problem (P) with the control constraints in (1.5) replaced by 

(1.7) ||it,(t)|| = l on [t,T-t] and ^ < ||u*(t)|| < ^ on [0,r) U (T - r,T], i = 1,...,m, 

where the time endpoint perturbation parameter r > 0 is arbitrarily small, and so (P’’) is not much 
different from (P). The purpose of the equality constraint relaxation on the small intervals adjacent to the 
time endpoints is to avoid degeneracy of necessary optimality conditions, which otherwise may hold for all 
the feasible solutions under some choice of nontrivial dual elements. Such a degeneracy phenomenon for 
necessary optimality conditions of the Pontryagin Maximum Principle (PMP) type has been discovered 
and well investigated in control theory with inequality state constraints; in particular, for Lipschitzian and 
compact-valued differential inclusions as in [4, 33, 42]. Our case is significantly different from the previous 
studies in both directions of the problem setting and the results obtained. We derive nondegenerate 
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necessary optimality conditions for intermediate local minimizers of {P'^) with t > 0 while the passage to 
the limit therein as r | 0 leads us to the conditions that generally degenerate. As examples show, even the 
degenerate optimality conditions obtained in this way for (P) can be useful to find optimal controls, but 
anyway we treat as our main result the more trustworthy ones established for (P’’). Of course, there is no 
difference between problems (P) and (P"^) if the equality constraints (1.5) are not imposed. 

Although the obtained necessary optimality conditions for the controlled sweeping process are con¬ 
structively expressed via the problem data and turn out to be efficient as illustrated by various examples 
presented in this paper, they are rather complicated and contain measures, which is not surprising for 
state-constrained systems. However, a crucial advantage of the method of discrete approximations and its 
strong convergence established below is that we can stop, with any prescribed accuracy, at a suitable step 
of discretization and treat the corresponding optimal solution to the discrete problem satisfying the (much 
simpler) discrete optimality conditions as an approximate/suboptimal solution to the continuous-time one. 

The rest of the paper is organized as follows. The main result of Section 2 shows that the class of 
absolutely continuous controls (u(-),6(-)) and the corresponding absolutely continuous trajectories a;(') is 
a right choice for feasible solutions to (P) and (P’"), since such a control pair satisfying the polyhedral 
constraints (1.4) ensures the existence of an absolutely continuous solution to the Cauchy problem (1.3) 
under an appropriate constraint qualification, which is also used in deriving necessary optimality conditions. 

The major aim of Section 3 is to construct well-posed discrete approximations of the optimal control 
problem (P’’) for any r G [0,T], with P° := P, such that they admit optimal solutions whose piecewise 
linear extensions on [0,T] converges to the given intermediate local minimizer (a;'^(-), u'^(-), fe'^(-)) of (P^) 
strongly in lT^’^[0,r] with some additional properties allowing us to derive nondegenerate necessary op¬ 
timality conditions for {■), u'^{■), b'^{■)) as r > 0 by passing to the limit from discrete approximations. 

This essentially distinguishes our new results in this direction from those obtained in the preceding paper 
[12] devoted to discrete approximations of the control sweeping process. As a crucial step of this procedure, 
we justify the strong IT^’^-approximation with the desired additional properties for any feasible solution 
to the controlled sweeping differential inclusion (1.3) without taking into account the cost functional (1.2). 

Since optimal control problems (P’’) for r G [0,T] and its discrete counterparts are intrinsically non¬ 
smooth due to the sweeping dynamics (1.3) and its finite-difference approximations, we need to employ 
suitable constructions of generalized differentiation satisfying extensive calculus rules to obtain necessary 
optimality conditions first for discrete-time and then for continuous-time systems. Section 4 is devoted 
to the description of such constructions and the explicit calculations of the major second-order one—the 
coderivative of the normal cone mapping—entirely in terms of the given data of the controlled sweeping 
process. These second-order calculations are certainly of their own interest while playing a crucial role in 
the efficient realization of our approach to deriving necessary optimality conditions. 

Section 5 presents necessary optimality conditions for the discrete-time optimal control problems appear¬ 
ing in the discrete approximation procedure for (P'^), r G [0,T], developed in Section 3. These conditions 
are obtained by reducing the discrete-time problems to nonsmooth mathematical programs with many ge¬ 
ometric and functional constraints with the usage of generalized differential calculus and the second-order 
calculations from Section 4. The conditions obtained are expressed explicitly via the problem data. 

Section 6 is a culmination of the paper. It contains the formulation and proof of the main nondegenerate 
necessary conditions for intermediate local minimizers of the sweeping control problem (P’’) whenever 
T G (0,T), their limiting versions as r 4- 0, and those for some special cases. The proof of the main result 
is rather involved and significantly depends on the major results obtained in the previous sections. 

The concluding Section 7 contains some applications to problems of quasistatic elastoplasticity with 
hardening and also present several examples showing the strength and illustrating specific features of the 
necessary optimality conditions obtained for the controlled sweeping process. 

The notation of this paper is standard in variational analysis and optimal control; see, e.g., [24, 42]. 
Recall that B stands for the closed unit ball of the space in question, B{x, r) := x-\-r'E, and N := {1, 2,...}. 
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2 Feasible Solutions to the Controlled Sweeping Process 

To begin our study, we want to make sure that the choice of absolutely continuous controls in problems {P'^) 
as T £ [0, T] is appropriate from the viewpoint of feasibility, i.e., such a choice of (m(-)i b(-)) in (1.4) ensures 
the existence of a solution x(-) to the Cauchy problem in (1.3), which is at least absolutely continuous 
on [0,T]. Observe that the unbounded polyhedral moving set C{t) generated by such a pair («(•),&(•)) 
in (1.4) is not absolutely continuous on [0,T] in the Hausdorff sense for set-valued mappings (not even 
talking about Lipschitz continuity), and hence we cannot deduce the existence of an absolutely continuous 
trajectory x{-) of (1.3) from known existence theorem for the sweeping process; see, e.g., [13]. 

Consider first the general sweeping process (1.3) generated by an arbitrary closed and convex moving 
set C{t) in K", which is assumed to be nonempty for all t £ [0,T]. Denote by v{t) := 7rc(t)(0) the unique 
projection of the origin onto C{t) and define the shifted set K{t) := C{t) — v{t). The following result has 
been recently proved in [12, Theorem 2.1] while being our starting point in this section. 

Lemma 2.1 (existence of absolutely continuous s-weeping trajectories for general moving 
sets). Let the projection v: [0,T] —>■ R" be absolutely continuous on [0,T]. We assume that for any 
positive numbers r,e there is a number S = S{r,e) >0 satisfying the estimate 

( 2 . 1 ) 

i—1 

for every collection of mutually disjoint subintervals 

i 

{[ai,Pi]\i=l,...j} of [0,T] with \Pi - ai\ < S. 

i^l 

Then there exists an absolutely continuous solution of the Cauchy problem in (1.3). 

We now use this result to establish an existence theorem for absolutely continuous trajectories of 
(1.4) generated by absolutely continuous controls in the polyhedral description (1.4) under the Linear 
Independence Constraint Qualification (LICQ). Note that this qualification condition was missed in the 
statement of [12, Corollary 2.2], where the proof was given only in the case of m = 1 in (1.4). We are 
very grateful to Alexander Tolstonogov for observing that an additional condition is needed for the validity 
of the latter result for m > 1 and that the case of m = 1 follows from his more general recent existence 
theorem in [41]. The following result new for the case of m > 1 is what we needed to justify well-posedness 
of the absolutely continuous framework for feasible solutions to problems {P'^) as r £ [0,T]. 

Theorem 2.2 (existence of absolutely continuous sweeping trajectories for controlled polyhe- 
dra). Let the controls (u(-) = (ui(-),...,Mm(')) b{-) = {bi{-),... ,bm{-)) be absolutely continuous on 

[0,T], let the inequality system in (1.4) be consistent, {i.e., C{f) ^ 0)) for all t £ [0,T], and let 

(2.2) the vectors {ui(t)| i £ Ix{t)'\ be linearly independent whenever x € C(t), t £ [0,T], 

where Ix{t) '.= {* £ {1,... ,rn} with {ui{f),x) = hi{t)'\. Then the corresponding Cauchy problem in (1.3) 
admits a unique absolutely continuous solution x(-) on [0,T]. 

Proof. Let us show first that the (well-defined) projection v{-) is absolutely continuous on [0,r]. Indeed, 
observe that for each t £ [0,T] the vector v{t) solves a positive-definite parametric program under the 
imposed LICQ (2.2). Then all the assumptions of Robinson’s stability theorem from [34, Theorems 2.1 
and 4.1] are satisfied. Taking into account that / appearing in the stability condition [34, formula (2.4)] is 
the Lagrangian of our quadratic program, x is our v{t), C is the first orthant, and the parameter p therein 
is our t, we conclude from the aforementioned stability condition that the modulus of continuity of v{-) is 
proportional to the modulus of continuity of the problem data with respect to t. This readily justifies the 
claimed absolute continuity of the projection function. 
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To deduce the existence theorem for the sweeping process under consideration from Lemma 2.1, it 
remains to verify the validity of condition (2.1) in this case. Recall that a multifunction T : [0,r] ^ M" is 
absolutely continuous if for any e > 0 there is some (5 (e) > 0 such that the implication 

lA -ai\<S (e) dff (r (oi), L (A)) < e 

holds for any finite collection of mutually disjoint intervals [oi, A] C [ 0 , T], where cIb refers to the Hausdorff 
distance. For r > 0 and t S [0,r] define Tr(t) := K{t) n rB. Then it follows from the definition of K{-) 
that Tr{t) = with the notation 

:= {Cj{t) — v{t)) n rB and Cj{t) := {a; S K"| {uj{t),x) < j = 1,... ,m. 

We claim now the validity of the estimate 

(2.3) dist ( 2 :,r(,'’)(t)) < r \\uj{t) - uAAII + lAO “ AAI + \{uj{t),v{t)) - {uj{s),v{s))\ 

for all s,t G [0, T], r > 0, j € {1,..., m} and z € rr'^^(s). Note that (2.3) is trivially satisfied for z € 
and hence we may assume that z ^ rr'^^(t). Now fix such j, z and show first that 

(2.4) dist (z, r^-^yt)) = 7 := {uj{t),z + v{t)) - bj{t) > 0. 

Indeed, (2.4) follows from the easily verifiable fact that z := z — "/Uj{t) is a unique minimizer of the program 

(2.5) mm{ ||z - yf \ y G Cj{t) - u(t)} 

by ||itj(t)|| = 1. Observing that v{t) G C{t) C Cj{t) and 7 > 0 gives us 7 < 2 {uj{t), z), and so 

ll^ll^ = ll^ll^ - 27 z)+'y^ = ||zf + 7 (7 - 2 z)) < ||zf < 

which means that z S rB and thus z solves not just (2.5) but also the program 

mm{||z-y||^| y GTi^\t)} 

over the smaller constraint set. This verifies (2.4). Further, it follows from z G rr'^^(s) that 

7 = {uj{t)-Ujis),z) + {uj{s),z) + {Uj{t),v{t)) -bj{t) 

< - Uj{s),z) + bj{s) - {Uj{s),v{s)) + {Uj{t),v{t)) - bj{t), 

which implies (2.3) by taking into account that ||z|| < r. Interchanging the roles of s and t in (2.3) gives us 

dff(rW(s),FW(t)) < r \\uj{t) - Uj{s)\\ + \b{t)-b{s)\ + \{uj{t),v{t)) - {uj{s),v{s))\ 

for the same s,t,r,j,z. On the other hand, the functions ruj{-), 6 (-), and {uj{-),v{-)) are absolutely 
continuous on [0,T] since Uj{-) and 6 (-) were assumed while v was shown to be such. This tells us that 
the multifunction Fr'^^(-) is absolutely continuous for any r > 0 and j G {I,..., m}. We want to derive 
from here that rr(-) is absolutely continuous for any r > 0 , which would follow from this property of the 
intersection mapping Fr^^) n Fr^^(-). Since this intersection is bounded, the desired fact follows from 

(2.6) (t) n int F),^^ (t) A 0 whenever t G [0, T] 

by [29, Proposition on p. 274]. To verify (2.6), observe first that the assumed LICQ and nonemptiness of 
C{t) yields the existence of a (time dependent) Slater point x{t) in the description of C'(t): 

x{t)) < bj{t) for all t G [0,T], j = 1,... ,m. 
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Recall the inclusion v(t) G Cj{t) for all j = 1,..., m on [0, T], which means that 

v{t)) < bj{t) whenever t G [0, T], j = 1,..., m. 


Define x^{t) := ux{t) + (1 — for t G [0,T] and v G [0,1] and deduce that {uj{t),Xv{t)) < bj{t) and 

S int rB as t £ [0,T], j = 1,2, v G (0, min{r, 1}). It yields Xv{t) G (t) fl int F^^^ (t) justifying (2.6). 

It is proved therefore that Fr(-) is absolutely continuous on [0,T] for any r > 0, which means that 
whenever r, e > 0 there is 6 {r, e) > 0 such that the implication 

lA - ail < ^ (d e) y]. dn (F^ (ai), F^ (A)) < e 

holds for any finite collection of mutually disjoint intervals C [0,T]. This verifies that 


i 


E 


max dist (z, K (Bi)) < 

zeK{ai)nrM 


max dist (z, iF (A) n rB) < y^ (F^ (a^) F^ (A)) < £ 

z&K(ai)nrM 

i—1 i—1 


and thus completes the existence part of the proof of the theorem. Uniqueness follows from a well known 
argument based on the convexity of the moving sets and Gronwall’s lemma. A 

The next example demonstrates that just the consistency condition C{t) A 0 is not sufficient for the 
existence of absolutely continuous trajectories in (1.3), (1.4) generated by C°° controls {u{-),b{-)) for m = 3. 


Example 2.3 (LICQ is essential for the existence of absolutely continuous sweeping trajecto¬ 
ries). Consider the controlled sweeping system (1.3), (1.4) in K.^ generated by the controls 

ui{t) := Cl, U2(t) := —ei, U3{t) ■= { — cost, — sint), bi{t) = 1 , &2(t) := — 1 , ^3(t) := — cost — sint 

on [ 0 , 7 r], which are obviously C°° functions on this intervals while LICQ (2.2) fails. Then we have 

= and .(i) = A- 

({I}x[l,oo) for 0 < t < TT [(lA) for 0 < t < tt. 

Due to the discontinuity of v(t), the assumptions of both Lemma 2.1 and Theorem 2.2 are not satisfied 
although C{t) A 0 for ah t G [ 0 , 7 r]. Observe further that the corresponding Cauchy problem (1.3) with 
x(0) = (1,0) cannot have absolutely continuous solutions, since even the requirement x{t) G C(t) on [0,T] 
is not met for continuous functions x{t) by the discontinuity of C(t). 

Remark 2.4 (Lipschitzian sweeping trajectories). It follows from [12, Proposition 2.3] that the 
Lipschitz continuity assumptions on v{-) and K{-) n rB in addition to those imposed in Lemma 2.1 ensure 
the existence of a unique Lipschitzian solution to (1.3) with an explicit estimate of its Lipschitz constant. 
Then the arguments in the proof of Theorem 2.2 allow us to conclude under these assumptions that any 
Lipschitzian control pair (u{-),b{-)) in (1.4) generates a unique solution to (1.3) with the aforementioned 
properties. Furthermore, if all the controls u{-) and b{-) in (1.4) are uniformly Lipschitzian on [0,T] with 
the moduli and Lf,, respectively, then we can prove the existence of a constant M > 0 dependent only on 
T, xq, Lu, and Lf, such that the sweeping process in (1.3) is equivalent to the bounded differential inclusion 

—x{t)GN{x{t)-,C{t)')C\M'E> a.e. tG[0,T] with a;(0) = ccq £ (7(0). 

Such a boundedness reduction was first established in [40] for the sweeping process with an absolutely 
continuous moving set C{t), which is not the case here, and then was further developed in [11] for (1.3), 
(1.4) with m = 1 under uniform Lipschitzian assumptions. In this paper, in contrast to [11], we prefer not 
to require the uniform Lipschitzness of feasible controls that led us in [11] to incomplete results on discrete 
approximations and necessary optimality conditions. 
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3 Well-Posed Discrete Approximations 


In this section we start developing a discretization approach to the study of the optimal control problems 
T G [0,r], which will hnally result in deriving robust necessary optimality conditions for their local 
minimizers. As the first step of this device, we construct here well-posed discrete approximations for an 
arbitrary feasible solution to (P^) and establish an appropriate strong convergence of such approximations. 

To begin with, let us represent the controlled sweeping differential inclusion (1.3) under the polyhedral 
constraints (1.4) in an equivalent differential inclusion form for extended trajectories. Consider the vectors 
u := (mi, ..., Urn) S M""*, b := (6i,..., bm) S M™, and z := (a;, u, 6) € M" x R""* x and define the 
set-valued mapping P: K" x R”™ x R'" R" by 

(3.1) F{z) :=—N{x]C(u,b)) with C(u, &) := {x € R"| (wi, x) < i = 1,..., m}. 

Then the sweeping differential inclusion in (1.3) under constraints (1.4) can be rewritten as 


(3.2) z(t) e G(z(f)) := P(z(f)) X R"™ X R"* a.e. t G [0,T], 

where the initial condition 2:(0) = (xq, rt(0), b{0)) is such that {ui{0), xq) < l>i(0) for alH = 1,..., m. We can 
see that form (3.2) treats controls b{-)) and the corresponding sweeping trajectories x(-) symmetrically. 
Theorem 2.2 tells us that any choice of absolutely continuous functions u(-) and b(-) satisfying LICQ (2.2) 
generates a feasible solution z(-) for the extended differential inclusion (3.2). 

Now the constrained sweeping system (1.3), (1.4) is written in the conventional form of the theory of 
differential inclusions with fixed right-hand sides as in [25, 42] while they do not possess major properties 
under which this theory has been developed. Indeed, the right-hand side of (3.2) is described by a highly 
irregular (discontinuous and unbounded) set-valued mapping. Further, it follows from definition (1.6) of 
the normal cone that (3.2) implicitly contains the state constraints on the trajectories z(t) given by 

(3.3) x{t) G C[u{t),b{t)) for all t G [0,P]. 


Moreover, besides the more conventional state constraints of the inequality type (1.4), each optimal control 
problem (P'^) as t G [0,r] contains those (nonsmooth and irregular) of the equality type given by (1.7), 
which has not been investigated in optimal control theory even in much simpler settings; see the discussion 
in Section 1. All of this emphasizes serious challenges we face to study these problems. 

The next theorem establishes the aforementioned strong -approximation of a given feasible solution 
to (3.2) subject to the state constraints in (1.7) (those in (3.3) are contained in (3.2)) by a sequence of 
feasible solutions to its discrete counterparts. The underlying difference of this theorem from the previous 
one in [12, Theorem 3.1] is deriving, under additional assumptions, new approximation properties that play 
a crucial role in the subsequent passage to the limit from optimality conditions for discrete-time control 
problems. We check then that the additional assumptions imposed are not actually restrictive. 

In what follows the symbol jrik) stands for the smallest natural number j such that tj>T whenever 
T G [0, r], while j'^{k) signifies the largest j G N with t^ <T — t. 

Theorem 3.1 (strong discrete approximation of feasible solutions). Let z(-) = (x(-),u(-),b(-)) G 
1F^’^[0,T] := 1F^’^([0,Tj;R”"*"”™’*'™) be an arbitrary feasible solution to problem {P'^) with any fixed pa¬ 
rameter T G [0,r] and define the uniform discrete partitions of [0,r] by setting 

(3.4) Afc := {O = tp < < ... < = P} with hk := tj_^^ — tj I 0, j = 0,... ,k — 1, as fc —>■ oo. 

Assume in addition that z{-) satisfies the following properties at the mesh points {all these properties are 
automatically satisfied if, e.g., z{-) G 1T^’°°[0,T]); the differential inclusion (3.2) holds for z{-) at all t^, 
j = 0,..., fc — 1 for each k gN, we have 


k-l 


1=0 


x(t)+i) - x{tk) 
kk 


-Htj) 

8 


2 

—I 0 as fc —>■ oo. 


(3.5) 







and there exists a constant M > 0 independent of k such that 


(3.6) 


k-l 

El 


- x{t]) ^ 


hk 


-m) 


< M, 


u{tl)-u{^) 


hk 


< M, 


6(t^) - 6(ig) 


hk 


< M, 




k-2 




3=0 


h{t]^2)-Wj+i) -h{t]) 


3=0 


< M 


whenever k G N. Then there is a sequence of piecewise linear functions z^{t) := {x'^(t), (t), (t)) on 

[0,T] with (a;^(0), M^(0), 6^'(0)) = (xq, m(0), 6(0)), 


(3.8) 


\\uHtj)\\ = 1 for j = jrik),... ,r{k), 

^ < ||Mf(ij)|| < ^ for j < jr{k) - 1 and j > f{k) + 1, i = l,...,m, 


satisfying for all k gN the extended finite-difference inclusions 

(3.9) x^{t) = x^{tj)(t — tj)vj whenever t^<t<ff_^i with € F(z^( 6^)) , j = 0, ...,fc —1, 
and such that the functions z^f) converge to z{-) in the norm topology of W^’^[0,T], i.e., 

T 

(3.10) z^{f)^z{f) uniformly on [0,T] and J \\z^{t) — z(t)\\^ dt ^ 0 as fc —?> cx). 


Moreover, there exists a constant M > M depending only on the total variation of uf) on [0,T] so that 
for every k gN we have the estimates 


(3.11) 




hk 


< M, 


b^(t\)-b'^{tl) 


hk 


< M, 


(3.12) var(u^; [0, T]) < M, and var(6^; [0, T]) < M. 

Finally, the sequence {x^(-)} has uniformly bounded variations on [0,T] being in addition uniformly Lips- 
chitzian on [0,T] if x{-) is Lipschitz continuous on this interval. 

Proof. Following the proof of [12, Theorem 3.1], for any fc € N we construct the continuous-time functions 
y^{t) := (yi(i), 1 / 2 ( 0 J 2/3 (0) the piecewise linear extensions on [0,r] of the discrete-time triples 

{ylit)),yl{t)),y'^[t))) := {x{t)),u{t^^),b{t^^)), j=0,...,k, 

and denote by w^{t) = (ref (t), w^it), w^{t)'^ := y^{t) their derivatives at non-mesh points. The assumptions 
in (3.7) yield vai{wf; [0, T]) < M for z = 2, 3 whenever k gN. It follows from the above that 

y^f) ^ z{-) uniformly on [0,T] and (•)—>• z(-) strongly in L^[0,r] as fc —> 00. 

Defining u^{t) := 2 / 2(0 'U^(O) = 2 / 2 ( 0 ) = "“(0) and deduce from (1.7) and the constructions 

above that the constraints on u^{-) in (3.8) and (3.11) hold. Fix k gN, denote tj := tj as j = 1, ..., k — 1, 
and construct the desired functions 6^(t),x^(t) on [0,T] by indnetion. To proceed, put (x^(0),6^(0)) = 
(xo, b{0)), suppose that the value of x^{tj) is known, and define 6f(t) at the mesh points so that 

b^dt3) ■■= {x\t,),u’l{t,))ii y^t,) = ( 2 /^( 0 ), 2 / 2 ^( 0 )), 

0^(0) > y^,{tj) > (2/^(0),2/2*(0))- 
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We can clearly arrange - h{tj) = “ 

Using the projection = TTF{x'‘{tj),u'^itj),b'‘{tj))iwij), define next x'^{t) on by (3.9), 

then construct b^iti^i) as above and extend it linearly to Observe that our construction yields 


bi (h+i ) - b- {tj) h {tj+i) - bi {tj ) 


b^i (^i+i) - b{tj+i) b\ {tj) - b{tj ) 

hk hk 


hk hk 


(3.13) 


X^{tj + i)-x{tj + i) 


x'^itj)-x{tj) 


< 


hk -^ / \ hk 

3 _ ryku _ ^f-h ^k( 


U^{tj 


X^{tjFl) — X^(tj) x{tjFi) — x{tj) X^{tj) — x{tj) 


for all j = 0,..., fc — 1. We also have the equalities 


(3.14) 


"ito+i) - x{tj+i) - x{tj) 


hk 


hk 


hk 


k - x{tj) 






= \\v^-w'^ 


ijW- 


Furthermore F{x^{tj),u^{tj),b'‘{tj)) = F{x{tj),u{tj),b{tj)), = 7r^(s(t^),a(t^.),5(t^))(wy), and 

_2 k _1 

“ ^ull = ;F{x{tj),u{tj),b(tj))^ 


(3.15) 


1=0 


3=0 

X{tj+i) -X{tj) 


3=0 


hk 


-x{tj) 


< M. 


Employing this together with (3.14) gives us 
(3.16) 


fe-i 

E 

3=0 


\\x'={tj+i) - x'^itj) 

xitj+i) - x{tj) 

hk 

hk 


< M, 


which readily implies the estimates 


(3.17) 


x'^itj)-x{tj) 


hk 




\x^{ti+i) - x{U+i) x^{ti)-x{ti) x^{to)-x{to) 


hk 


hk 


< M 


for every j = 0,..., fc — 1 as fc € N. As a result, it follows from (3.13), (3.14), (3.15), and (3.17) that 


E| 

3=0 


hk 


hk 


Combining the latter with (3.7), we arrive at 


k-2 

El 

2=0 




hk 


hk 


< M + Mva,r(u; [0, T]). 


< 3M + 2Mvar(u;[0,r]), 


with verifies the validity of (3.12) for b^. Observe simultaneously the fulfillment of the estimate for b^ in 
(3.11), which follows from (3.17) and the representation b^{ti) — bi{ti) = (x^{ti) — yi{ti),Ui(ti)). 

Next we justify the lU^’^-convergence in (3.10) for which it suffices in fact to check the L^-convergence 
of x^ and h^. To verify the former one, observe by (3.5) that 


pT k-l k-l -,,k \ _ ^ifk\ 

/ \\x\t)-w’i{t)rdt=Y,hk\\v’;-w>;r<Y,hk 

JO .n .-n 


3=0 j=0 

as k ^ oo. The claimed convergence for 5^ follows from (3.13) and 


^ 0 


r wb'^it) 

Jo 


k-l 


- Wo 


< 2 / \\x^it) -'w'l{t)\\‘^dt+ 2M‘^^hk\\ui{tj+i) - Ui(tj)\\" 

JO j=o 
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by the absolute continuity of u(-) on [0, T], The last statement of the theorem on {a;^(-)} follows immediately 
from (3.6), (3.16), and the fact that x(-) has bounded variation on [0,T]. To complete the proof of the 
theorem, it remains to observe that the validity of all the assumptions in (3.5)-(3.7) for the case of 
z(-) € 1T^’°°[0,T]) is a direct consequence of the definitions. A 

It is not hard to check that all the assumptions of Theorem 3.1 are satisfied if z(-) is piecewise C^[0,T] 
with z(-) € iiy([0,T]). In this case, more general than z{-) € 1T^’°°[0,T], the derivatives appearing in 
(3.2) and in the subsequent formulas are the right derivatives. 

Our next goal is to construct a well-posed discrete approximation for a given local optimal solution of 
the control problem {P'^) as r G [0,T'], which satisfies the assumptions of Theorem 3.1. We consider a 
rather broad class of local minimizers introduced and first studied for differential inclusions in [23] under 
the name of intermediate local minimizers (i.l.m.). This notion obviously covers strong local minimizers 
(corresponding to a = 0 in the definition below) and occupies an intermediate position between weak and 
strong minimizers in dynamic optimization and optimal control; see [23] and [25, Chapter 6] for more 
details. Note that a related class of local minimizers in optimal control problems for differential inclusions 
has been studied later under the name of W^’^-minimizers; see, e.g., [42]. We now present an adaptation 
of the i.l.m. notion to the case of the sweeping control problems (P’’) under consideration. 

Definition 3.2 (intermediate local minimizers for the controlled s-weeping process). Fix any 

T G [0,T] and consider a feasible solution z(-) = (a;(-),'u(-), 6(-)) G W^’^[0,T] to (P^). We say that z{-) 
is an INTERMEDIATE LOCAL MINIMIZER for this problem if there are numbers a > 0 and e > 0 such that 
J{F\ < J[z\ for any feasible solution z{-) = (a;(-), u(-), 6(-)) to {P'^) satisfying 

II (x(t), u{t),b{t)) - {x{t), u{t),b{t)) II < e as t G [0, T], and 


(3.18) 




i{t) 



b{t) - b{t) 


dt < e. 


It is easy to see that the general setting of a > 0 in (3.18) reduces to the cases when either a = 1 or 
a = 0, and that from the viewpoint of necessary optimality conditions it suffices to examine only the case 
of intermediate local minimizers with a = 1, which we do in what follows. 

Given a number r G [0, T] and an i.l.m. z(-) for (P^), construct now the family of discrete approximation 
problems (P^), fc G N, having optimal solutions that converge to z(-) in some strong sense, which eventually 
allows us to derive eventually necessary optimality conditions for z{-) by passing to the limit from those 
in discrete approximations. Suppose that for the given i.l.m. z(-) and the discrete mesh in (3.4) all the 
assumptions (and hence conclusions) of Theorem 3.1 are satisfied and then define each problem (P^) by: 


fc-i 


minimize Jk[z^]:= g}{x%) + hk^^J {t ), x ], u'f , b ’], 


y,k ^,k ^,k ^k _ 


(3.19) 


-f 


1=0 

fe_i q+i , 


hk 


hk 


hk 


AC —1 -J-hi / 

E / 

i=n xjf ^ 


1=0 

■ dist^ ^ 


A+i 


hk 


-x{t) 


+ 


^1-1-1 “l 


hk 

k-2 


-u{t) 


2 

A+1 A 


hk 


b{t) ^ dt 


hk 


(—oo, M]J + dist^ ^ 


1=0 


‘‘ ''»|,(-oo,M|) 


hk 


+ dist^ ^ 

1=0 

k-2 

-l-dist^ ^ ^ 


fe-2 ,, ,.k 


'»l-^2 ^'^1-H A ^1 


■ n ^k 

J=0 

^“2 ik _ nuk Lfc 

Oj+2 ^ 0 ^ + 1 + ^1 


1=0 


hk 


{-oo,M]^ 
, (-oo,M]^ 
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over collections := (xq, ..., x^, Uq, ..., bg,..., b\) with := ..., G K”"* for every index 

j = 0,... ,k subject to the constraints in (3.8) together with 

(3.20) G Xj + hkF{x’^,u’^,b'^) for j = 0,..., fc - 1 with (xQ,Uo,bQ) = (xo,u(0),b(0)), 

(3.21) (uki,Xi/,} < b^i for i = 

(3.22) II (x^, b^) - (x(t^), u(t^), b(t^)) II < e/2 for j = 0,...,k, 


(3.23) 


(3.24) 







hk 



hk 



hk 


< M + 1 , 


b>l-b^ 

hk 


<M + 1, 


(3.25) 


fe -2 


E 


u ’;+2 - + u ) 


k-2 

< M + 1, and 

3=0 


b)^,-2b]^,+b] 


< M + 1, 


where M is taken from in Theorem 3.1 while e > 0 is taken from Definition 3.2 with 0 = 1 . Note that the 
index j plays the role of the discrete time in (P^) and that inclusions (3.20) correspond to those in (3.9) 
at the mesh points of A^. Observe that, in contrast to part (3.8) of the state constraints, the other part 
(3.21) needs to be imposed only at the endpoints {x/.,u/.,b//) while the counterparts of (3.21) at {xj,Uj, bj) 
for j = 0,..., fc — 1 follows from (3.20) due to the structure of F in (3.1) by the definition of the normal 
cone in ( 1 . 6 ). It is important to emphasize that the set of feasible solutions to each problem (P^) with 
T G [OiT] and fc G N sufficiently large is nonempty by Theorem 3.1. 

To employ and justify the method of discrete approximations in deriving necessary optimality conditions 
for the control sweeping process, we need to make sure that for all r G [0, T] and all fc G N sufficiently large 
each problem ) admits an optimal solution. Despite the finite-dimensionality, this issue is nontrivial for 
{Pf) due to the possible nonclosedness of the feasible solution set to this problem because of the dynamic 
constraints (3.20) generated by the normal cone to the moving set in (3.1); see [12, Example 4.5]. To 
overcome such a possibility, we employed in [12] the Positive Linear Independence Constraint Qualification 
(PLICQ) for the given i.l.m. z{-) in the original problem (P) formulated as follows: 


(3.26) 


OilUlit) = 0, tti > 0 

{x{t) ,u{t) 


0 ^ = 0 as i G l[x{t),u{t),b{t)) on [0 ,r], 


where the collection of the active constraint indices i G I{x{t),u{t),b{t)) is defined by 


(3.27) 


I{x,u,b) := {z G {1,... ,m}| {u^,x) =5*}. 


This condition, which permits the linearly dependence of active gradients, is obviously weaker than the 
classical LICQ for z(-) on [0,T] that has already be used in Theorem 2.2. It is worth mentioning that in 
our polyhedral setting (1.4) under the additional assumption of C{t) 7 ^ 0 it follows that (3.26) corresponds 
to Slater’s condition, while we keep the term PLICQ here since a more general framework of moving sets 
described by smooth inequalities will be considered in Section 4 and further developments. 

In [12, Theorem 4.4] we proved the existence of optimal solutions to problem {Pk) with Pk := P^ for all 
large numbers fc G N under the validity of PLICQ (3.26) by using the normal cone/subdifferential structure 
of the constraints in (3.20) and employing Attouch’s theorem on the subdifferential convergence (see, e.g., 
[ 6 , Theorem 12.35]) as well as the extremal principle of variational analysis from [25, Theorem 2.8]. The 
given proof holds without any change in the case of problems {Pjf ) with r > 0 by taking into account the 
existence of feasible solutions to (Pjf) justified in Theorem 3.1. This brings us to the following result. 


12 



































Proposition 3.3 (existence of discrete optimal solutions). Let the cost functions (p and i{t,-,-) be 

lower semicontinuous around the given i.l.m. z{-) satisfying PLICQ (3.26) whenever t £ [0,T]. Then for 
all T £ [0,T] and all k sufficiently large there exist optimal solutions to the discrete problems (P^)- 


To proceed next with establishing the concluding result of this section on the desired strong convergence 
of optimal solutions for {P^) to the given local minimizer z{-)^ we need to impose one more requirement 
on z(-). Fix any quadruple (t, x, u, b) and denote ipif, x, u, b, u, w, v) the convexification of the integrand in 
(1.2) on the set F(x^u,b) from (3.1) with respect to the relocity variables {v,w^v), i.e., the largest convex 
and lower semicontinuous (l.s.c.) function majorized by £{t, x, u, b, •, •, •) on this set. Then for all r > define 
the relaxation {K^) of problem {P'^) as follows: 


(3.28) 


minimize 


J[z]-.= ip{x(T)) + f iF{t,x{t),u{t),b{t),x{t),u{t),b{t)) dt 

Jo 


over the triples z{t) = {x{t), u{t), b{t)) of absolutely continuous functions on [0, T] satisfying the constraints 
in (1.3)-(1.5) and (1.7). It follows from the construction of ip and the convexity of the set on the right- 
hand side of (1.3) that the relaxed problem (i?^) reduces to the original one (P^) if the integrand i in (1.2) 
is convex and l.s.c. with respect to the velocity variables {v, w, v). In the general case we say that z(-) is a 
relaxed intermediate local minimizer (r.i.l.m.) for (P^) if it is an i.l.m. for this problem with J[z\ = J[z\. 

A remarkable phenomenon well-recognized in control theory for continuous-time systems reveals that 
in many nonconvex settings the value of the cost functional does not change under the integrand convexi¬ 
fication with respect to velocity variables. It is known as “hidden convexity” being related to Bogolyubov- 
type relaxation results and Lyapunov’s convexity theorem for integrals of set-valued mappings; see, e.g., 
[7, 25, 42]. To the best of our knowledge, the most general Bogolyubov-type theorem is obtained in [14] 
for optimal control problems governed by differential inclusions satisfying the so-called “modified one-sided 
Lipschitzian” condition with respect to state variables. However, the latter condition does not hold for 
the sweeping inclusion (1.3). Thus we cannot drop so far the relaxation property of intermediate local 
minimizers in the following theorem, which is crucial for justifying the method of discrete approximations 
to derive necessary optimality conditions for the sweeping control problems under consideration. 


Theorem 3.4 (strong convergence of discrete optimal solutions). Given an arbitrary number 
T £ [0, T], let z{-) = {x{-), u(-), b(-)) be a r.i.l.m. for problem (P”) satisfying the assumptions of Theorem 3.1 
and Proposition 3.3 and such that tp is continuous around x{T), i{t,-,-) is continuous around {z(t), z{t)) 
uniformly on [0,r] while £(-,z,z) is a.e. continuous on [0,r] being uniformly majorized by a summable 
function near the given local minimizer. Then any sequence of piecewise linear extensions of the optimal 
solutions z^ = {x^,u^,b^) to the discrete problems (P^) converges to z{-) in the norm topology of W^’^[0,T]. 
Furthermore, we have the following estimates: 


(3.29) 


lim sup 

k-f-oo 


u\ - Uq 


hk 


< M, 


lim sup 

k—foc 



< M for all k gN, 


fc-2 ,, „-,k 


(3.30) 


lim sup 

k—¥CO r, ' 


wj+2 - + u] 


hk 


<M, and limsup^ — ^ ^ 


k—^oc 


3=0 


hk 


< M, 


where the constant M > M is taken from (3.12). 

Proof. Fix any sequence of the (well-defined by Proposition 3.3) extended optimal solutions z^{-) to 
problems (P^) and observe that strong LF^’^-convergence to z{-) on [0,T] as well as the properties (3.29) 
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and (3.30) follow directly from the equality 


/ r'~^ 2 2 2 

Jim ij ^ x^{t) — x(t) + u^{t)—u{t) + b^(t) —b{t) ^ 


(3.31) 


- dist^ ^ 


k-2 


dt 


, (- 00 , M]^ + dist^ ( X! I ,(-oo,^]) 


i=o 


hk 


—2Ti^ +7“/^ 


i=o 


hk 


_. — 2 /)'? + ft'? _\' 

, (-oo,M]^+dist^(^^ -_2±i-1 ,(-cx),M]^ 


i=o 


= 0 


due to the initial conditions (a;''(0), u''(0), ft''(0)) = (a;(0), m( 0), &(0)) as fc G N. To justify (3.31), suppose 
the contrary, i.e., the limit along a subsequence therein (no relabeling) equals to some 7 > 0. By the weak 
compactness of the unit ball in L^[0,T] := i^([0,T];R"' x R"™ x R™), find (u(-), w(-), z/(-)) G T^[0,T] and 
(if necessary) another subsequence of {z’^{-)} so that 

(x'‘(-),u^(-),P(-)) (v(-),w(-),iy(-)) weakly in L^[Q,T]. 

Next we define the absolutely continuous triple z{-) := (x(*), 6 (-)): [0,T] —)■ by 


t 

z{t) := (a;o,M(0),6(0)) + j {v{s),w{s),iy{s)) ds, 

0 


tG [ 0 ,T], 


for which z{t) = {v{t),w{t),i'{t)) a.e. on [0, T], Applying Mazur’s weak closure theorem gives us a sequence 
of convex combinations of {x^{-), u^{-), converging to (r’('), ^^(■)! ^(O) strongly in L^[ 0 ,T] and thus a.e. 
on [0, T] along a subsequence. Then passing to the limit as fe —i> oo in the discrete inclusions (3.20) with the 
convex normal structure of the mapping F from (3.1) and employing the aforementioned Attouch’s theorem 
tell us that x{-) satisfies the sweeping inclusion (1.3) with the set C{t) in (1.4) defined via u{-) and ft(-). 
The validity of the r-constraints in (1.7) for u(-) follows, whenever r G [0, T], from the uniform convergence 
on [0,T] of the designated sequence of convex combinations of u^(-) to the limiting control function m(-). It 
also follows from the strong L^-convergence of the above convex combinations of {x’^{-), b^{‘)) that the 

limiting triple z(-) belongs to the prescribed e-neighborhood (in W^'^) of the i.l.m. z{-) from Definition 3.2. 

It remains to pass to the limit in the discrete cost functional (3.19) along the optimal triple z^{-) for 
{PI) as k ^ oo. We can directly deduce from the construction of ip and its convexity in velocities that 


fc-i 


ipit, x{t), u{t), b{t), x{t), u{t),b{t)) dt < hmmf hk , 

' ' fc—700 z / 

3=0 




u ft 


p.k 

Xj + 1 


-x] 


u'j+l 




hk 


hk 



By the structure of (3.19), the lower semicontinuity of the total variation, and the choice of 7 above we get 


(3.32) 


J[z\ -I- 7 = ip[x{T)) + / iF(t,x{t), u{t),b{t),x{t),u{t),b{t)) dt + ^< liminf Jk[z 


k—¥oo 


by using the Lebesgue dominated convergence theorem due to the assumptions made. On the other hand, 
applying Theorem 3.1 to the local minimizer z(-) under consideration gives us a sequence {^^(•)} of the 
feasible solutions to {Pk) that approximates z(-) in the norm topology of 1T^’^[0,T]. Since Zk{-) is an 
optimal solution to problem (P^) while z'^{-) is feasible to it for each k, we have 


(3.33) Jk[z^] < Jkiz’^] whenever fc G N. 

It now follows from the structure of the cost functional (3.19) in (P^) with M > M, the strong 
convergence in Theorem 3.1, and the assumed continuity of ip and ^ that Jk[z^\ J[z\ as /c —)• 00 . Thus 
by taking (3.33) into account we obtain 


(3.34) 


limsup Jfc< J[z\. 

k—¥co 
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The obtained relationships (3.32) and (3.34) together with the assumption on 7 > 0 imply that J\z\ < J\z] 
contradicting therefore the choice of z(-) as a r.i.l.m. for {P'^). Hence 7 = 0, which shows that (3.31) holds 
and thus completes the proof of the theorem. A 

4 Generalized Differentiation and Second-Order Calculations 

After establishing well-posedness of the discrete approximation problems (P^) and the desired strong 
convergence of their optimal solutions to the given r.i.l.m. z(-) for the sweeping control problem (P^) 
with any fixed r G [0,r], our further strategy is as follows: obtain necessary optimality conditions for 
finite-dimensional discrete-time problems (P^ ) whenever fc € N and then justify the possibility of passing 
to the limit as A: —>■ 00 in the obtained discrete relationships as to derive necessary optimality conditions 
for z(-) in (P’’). Since problems (P^) and (P’’) are always nonsmooth due to the dynamic constraints 
independently on the smoothness of the cost functions ip and i in (1.2), we have to employ appropriate 
generalized differential constructions of variational analysis enjoying comprehensive calculus and robustness 
properties. In our setting not only first-order but also second-order generalized differentiation is needed. 

The main results of this section give upper estimates as well as precise formulas for calculating the 
coderivative of the normal cone mapping to moving sets as in (3.1), which is a second-order object playing 
a decisive role in the subsequent results of this paper. We begin with some basic definitions from generalized 
differentiation while referring the reader to [24, 36] for more details on the first-order constructions and to 
[24] and the papers mentioned below for the second-order ones and their equivalent descriptions. 

Recall that, for a set-valued mapping/multifunction F: R" R"*, the symbol 

(4.1) LimsupP(a;) := {y G R"*] 3 sequences Xk —> I, yk —> y with yk G F{xk) for all fc G N} 

X—¥X 

signifies the (Kuratowski-Painleve) outer limit of F at x. Given a subset H C R” locally closed around 
a; G H, the normal cone to H at a; (known also as the limiting/basic/Mordukhovich one) is defined by 

(4.2) N{x;ft) = Nn{x) := Limsup {cone[a; — n(a;; fl)] } 

X—^X 

via the outer limit (4.1), where n(a;; fl) stands for the Euclidean projection of x onto H, and where ‘cone’ 
denotes the conic hull of the set. When H is convex, the normal cone (4.2) reduces to the classical one 
of convex analysis, while in general the cone (4.2) is nonconvex even for simple sets fl, e.g., for := 
{{xi,X 2 ) G R^l X 2 = |a;i|}. Nevertheless, the normal cone and associated subdifferential and coderivative 
constructions for functions and multifunctions enjoy full calculi based on variational principles; see [24, 36]. 
Given a set-valued mapping F: R" x R"* whose graph 

gphE := {{x,y) G R" x R™| y G F{x)} 

is locally closed around (x,y), the coderivative of F at {x,y) is defined by 

(4.3) i3*F(x,y)(u) := {uGR"| (u,-u) GiV((S,y);gphE)}, u G R™, 

where y = F{x) is omitted if F is single-valued. When F: R” R"* is smooth around x, we have 

D*F{x){u) = {VF{xyu} for all u G R", 

with A* standing for adjoint operator/matrix transposition of the Jacobian A = VE(a;). 

For a l.s.c. extended-real-valued function p : R” —> R with the domain and epigraph 

dom(^ := {a; G R"| p{x) < 00} and epiy? := {{x,y) G R"'''^| /i > ^^(a;)} 

its (first-order) subdifferential at a; G dom(/3 is generated by (4.2) as 

(4.4) dp{x) := {u G R"*! {v, —1) G IV((a;, p{x)-, epit^)}. 
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Our main objects here are the second-order generalized differential constructions defined by the scheme 
of [24] as follows. Given v S dip{x) from (4.4), the second-order suhdifferential (or generalized Hessian) of 
tp aA, X relative to v is the mapping d^(f{x, v): R” R" with the values 

(4.5) d^(p{x,v){u) := {D*d(p){x,v){u), u e R". 

Having an extended real-valued function p-. R" x R"^ —R of two variables (x, w) € M" x R'^ and its partial 
(in x) first-order subdifferential mapping 

dxp{x,w) := {set of subgradients v of pw '■= at x}, (x,ty) G domp, 

define the partial second-order subdifferential of in x at (Xjtu) relative to x £ dxp{x,w) by 

(4.6) d‘lp{x,w,v){u) ■.= {D*dxp){x,w,v){u), m G R". 

Note that for C^-smooth functions p the constructions in (4.5) and (4.6) reduce, respectively, to 
d‘^p{x){u) = {Vl^p{x)u}, d^p{x, w)(m) = { {Vl^p(x, w), Vl^p(x, w))}, uG R", 

expressed in terms of the classical (symmetric) Hessian matrices. The partial second-order construction 

(4.6) has been studied in [27] under the name of “extended partial second-order subdifferential” with 
the notation d^p. Since no other partial second-order constructions are used here, we drop both the 
word “extended” and the tilde-notation for (4.6). Our goal is to estimate and ealculate this second-order 
construction for the special class of functions arising in the controlled sweeping process (1.3). 

To proceed further, consider the smooth parametric inequality system 

(4.7) S{w):={xGM.'^\g{x,w)GM.^}, 

where R™ is the nonpositive orthant of R"* and g: R” x R*^ —>■ R™ is an arbitrary C^-smooth vector function. 
Associate with (4.7) the normal cone mapping Af: R” x R'^ R" defined by 

(4.8) Af{x,w) := N[x; S{w)) for x G S{w) 

via the normal cone (4.2) to the moving set S{w)^ which we denote as Ns(w){x) for convenience. It is easy 
to see that the mapping Af in (4.8) admits the composite representation 

Af{x,w) = dxp{x,w) with p{x,w) := (Jr™ og^(x,w) 

by using the indicator function ^Rm of the orthant R™. Thus we get by definition (4.6) that 

(4.9) d^p{x,w,v)(u) = D*Af{x,w,v)[u) for any vGAf{x,w) and u £ R". 

Dealing with the moving set (4.7), we use in what follows the coderivative form (4.9) of the partial second- 
order subdifferential of the function p in question. 

The main issue is to evaluate this construction entirely via the given data of (4.7). The next lemma 
based on the second-order chain rules from [26, 27] plays an important role in the subsequent calculations. 
Recall that the mapping M: R'* R** is calm at (s, q) G gph M if there are numbers /i > 0 and 77 > 0 with 

(4.10) M(s) n (q -I- 77B) C M(s) -I- 7i|[s — s||B whenever s G s-\-pM. 

Lemma 4.1 (coderivative of the normal cone mapping for smooth inequality systems). Let 

(x, w) G R" X R'^ be such that g(x, w) G R™, let v £ Af(x, w), and let 

I{x,w) := {i G {l,...,m}| gi{x,w) = 0} 

be the collection of active indices for (4.7) at (x,w). The following assertions hold: 
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(i) Assume that the partial gradients {Vxgiix,w)\ i G I{x,w)} are positively linearly independent and 
that the mapping d {{x,w,p)\ {g(x,w),p) + 1 ) G gphiVRm} is calm at (0,x,w,p) for all p G N-Rm(^g(^x,w)) 
with yxgix^wYp = X. Then for all u G R" we have the upper estimate 


D*J\f{x,w,v){u) C 

u { 

pGNM,m{g{x,w)) L 
^xg{x,w)*p—v 

(ii) Assume that the partial gradients {Vxgiix,w)\ i G I{x,w)} are linearly independent, and let the 
vector p G R™ he uniquely defined by 

p G Ns:p[g{x,w)), Vxg{x,w)*p = v. 


XX 

V72 


{p,g){x,w) 
{p,g){x,w) 


u + Vg{x, w)*D*NRrn (^g{x, w),p) {Vxg[x, w)u) 


Then for all u G R" we have the precise coderivative formula 


D*J\f{x, w, v){u) 


VL(P)5)(a;,w) 

'^lw{P,g)(.x,w) 


u + Vg{x, w)*D*Ns.rn (^g{x, w),p) {Vxg{x, w)u). 


Proof. We derive (i) from the second-order chain rule of the inclusion type established in [26, Corol¬ 
lary 3.2(b)], where the first-order qualification condition follows from the positive linear independence of 
the active constraint gradients in (4.7). Assertion (ii) is a direct consequence of the precise (equality type) 
second-order chain rule obtained in [27, Theorem 3.1] under the full rank condition, which is ensured here 
by the assumed linear independence of the active constraint gradients. A 

Next we apply these results to the case of bilinear vector function g(x,w) in (4.7), which covers our 
controlled sweeping setting in (1.3). Define 


(4.11) g{x,w) := Ax — b with w := {A,b) for x G R" and 6 G M™, 


where A is an m x n-matrix, and both A and b are variable. In this case system (4.7) is written as 
(4.12) 5'(A,6 ):={xGR”| Ax<6}. 

Taking into account that the values of S{-, ■) are polyhedral sets, we refer to (4.12) as to the polyhedral 
system. Note that the graph of S may not be a convex polyhedron in R" x R"™ x R™. For any fixed (A, b) 
the active index set from Lemma 4.1 reduces to 


I{x,A,b) := [i G {1,... ,m}| A^x = &i}, 

and we label {A^l i G I{x, A, b)} as active rows. Based on Lemma 4.1 and the affine structure of (4.11), we 
arrive at the next lemma, which relates the coderivative of Af with that of Agm. 

Lemma 4.2 (coderivative of the normal cone mapping for polyhedral systems). Let (x,A,b) G 
R" xR"™ xR™ he such that Ax < b, and let v G A/"(x, A, b) for the corresponding normal cone mapping (4.8) 
generated by the polyhedral system (4.12). Assume that the active rows {Ai\ i G /(x,A, 6)} are positively 
linearly independent. Then we have the upper estimate 


D*M{x, A,b,v){u) C 




( A*q 

\ 


' 



Piu + qix 




u 





p G iVRm(Ax — 6), A*p = V, q G D*Ns_m(^Ax — b,p)(Au) 



PmU + qraX 





< 

\ 

) 


> 
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for all u G R". If moreover the active rows {Ai\ i S I(x,A,b)} are linearly independent, then we have the 
precise formula for the coderivative calculation 


D*Af{x, A, b, v){u) 



✓ 

( ] 


' 

u< 


PiU + qix 


q G D* NB.yn{Ax — b,p){Au) 



PmU + qmX 

\ y 




weK”, 


where the vector p G is uniquely determined by A*p = v. 


Proof. Applying Lemma 4.1 to g{x,w) from (4.11) and using x"^ for the corresponding vector row yield 

Vg {x,A,b) = 


/ 

-T 

X 

0 

0 

\ 

A 

0 


0 


V 

0 

0 

-T 

) 


^lx{P^9)=0, ^l,(A,b){P,9) = {Pil \ \PmI\0)* 

for any fixed {x, w) with w = {A, b) and any p G R™. Observe that the mapping 

M(d) ■= {(x, w,p)| {Ax — b,p) + d e gphiVRm}, d = {di,d2) G R^"*, 

is automatically calm at {0,x, A,b,p) for any p as required in Lemma 4.1. This is a consequence of 
the polyhedrality of M by the classical Robinson theorem from [35]. Thus the asserted formulas follow 
immediately from Lemma 4.1 and the Jacobian and Hessian calculations given above. A 

Now we are ready to derive from Lemma 4.2 the desired results for evaluating the coderivative D*M of 
(4.8) entirely via the given data of (4.12) by using the calculations of D*N’gm available in the literature. 
Consider the mapping T{x, A, b) := —Af{x, A, b), which actually appears in the sweeping inclusion. 


Theorem 4.3 (coderivative of the normal cone mapping via the given data). In the setting of 
Lemma 4.2, suppose that the active rows {Ai\ i G I{x,A,b)} are positively linearly independent. For all 
u G R" and p G R"*, respectively, define the sets 

P{u) := {p G NRm{Ax — b) I A*p = —u} if u € Aj' and P{u) := 0 otherwise, 

{i\pi>0} 

Qi =0 if AiX < bi or if AiX = h, pi = 0, AiU >01 
qi>Q if AiX = bi, Pi = 0, AiU <0 / ' 



Then for all u G R" we have the upper estimate 


(4.13) 


D*T{x,A,b,v){u) C IJ < 


p&P{u) 

qeQ{p) 


' 

( \ 

' 


qix-piu 



qmX - PmU 



\ -9 y 

y 


. 


If furthermore the active rows {Ai\ i G I{x,A,b)} are linearly independent, then either 



f 

( ^ 




qix-piu 


(4.14) D*P{x,A,h,v){u) = y < 




qGQip) 


qmX - PmU 




\ -9 y 

y 


if ue Pi [Ai\ , 

{i|pi>0} 


or D*F{x, A,b,v){u) = 0 otherwise. Here the vector p G Ns.'^^Ax-b) uniquely defined by A*p = —v. 
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Proof. Observe that D*T{x, b, v)(u) = D*JV{x, A, b, —v){—u). The claimed results follow from Lemma 4.2 
by substituting therein the precise coderivative calculation 


£)* (a,/3)(7) = 


if /3i7i ^ 0 for some i, 


{rj € IR™| r]i = 0 Hi € Ii and rji > 0 for i € J2} otherwise 
given in [16, p. 1215], where the index subsets of {1,..., m} are defined by 

Ii := {f| Oi < 0} U {i| ai= /3i = 0, 7^ < O}, I 2 ■= {i| = Pi = 0, 7* > O}. 

This verifies both coderivative formulas (4.13) and (4.14) of the theorem. 


A 


5 Necessary Optimality Conditions for Discrete Approximations 

The aim of this section is to obtain necessary conditions for optimal solutions of the discrete approximation 
problems [P^) for any fixed r G [0,T] and fc G N. First we derive optimality conditions for a generalized 
version of (P^ ), where the dynamic constraints (3.20) are described by an arbitrary closed-graph mapping 
F. Then, by using the coderivative calculations of Section 4, we arrive at optimality conditions expressed 
entirely via the problem data of [P^) with F given in the particular normal cone form (3.1) of the sweeping 
process under consideration. Our standing assumptions in this and next sections are that the cost functions 
ip and i(t, •, •, •, •, •) are locally Lipschitzian around the points in question. Note that the subdifferential (4.4) 
of the running cost £ is taken with respect to its all but t variables. In what follows we drop indicating the 
time-dependence of £ for brevity and use the notation 

(5.1) [c, q] := (ciQi,Cmqm) G R”™ and rep„(x) := (x,..., x) G M”™ 

for vectors c = (ci,..., Cm) G M*", x G R", and g = (gi,..., qm) G R"™ with gi G R" as i = 1,..., m. 


Theorem 5.1 (necessary optimality conditions for general discrete inclusions). For fixed fc G N 
and r G [0, T] let = (xo, x^ ..., x^, Hq, ..., 6g,..., &^) he an optimal solution to the discrete problem 

(3.19)-(3.25) written in the format of {Pjf) but with the discrete inelusion (3.20) governed by a general 
elosed-graph mapping F. For each j = ,k — 1 we denote 


(5.2) 





dt. 


Then there exist dual elements > 0, G R™, G pj = G 

■^n+nm+m j = Q,..., fe and subgradieut vectors 


(5.3) 




^3 ’ ^3 


,vf,vf) €d£ 



j = 0,... ,fc - 1, 


such that the following conditions are satisfied: 


(5.4) 


k-1 

A'= + ||a1l + ni+Ellpfll 


3=0 


boll + boll 


(5.5) 


a^{{ul,,xt)-bl)=0, 


i = 1,..., m, 


(5.6) G fV(||Uji||; [1/2, 3/2]) /or j = 0,.. .,j^(fc) - 1 and j = f (k) + I,..., k, f = l,...,m, 
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(5.7) -pfeX>^d^{xl)+J2a>^ul,-, pf = -[a\Tep^{xt)]-2[e„4], pf = a\ 


2 = 1 


(5.8) 


(5.9) 


pT+1 = + K Pj+i = ^ )’ J = 0,..., ft - 1 


C xk xk uk uk 

Pj + l Pj ^k^..xk Pj+1 Pj 

~~h,. ^ -—h, ^ 


bk _ bk 

k^,,uk ^i+1 xj \k^,,bk ^xk \k 

’ Jlk ^ ’^i+1 ~ ^ 




e (o, il [«‘, a*], 0, o) + w ((i‘, a*, if, ; gph f) , j = o,..., a - i. 


hk 


Proof. Throughout the proof we omit indicating the (fixed) upper index ‘ft’ from the statement of this 
theorem; the dependence of the result on ft will be needed in Section 6. Let 

y '■= (a^Oi ■■■ ,Xk,uo, ■■■ ,Uk, bo,... ,bk,Xo, ■■■, Xk-i,Uo,... ,Uk-i, Bq, . ■., Bfe-i), 

where xq is fixed. Take £ > 0 from the construction of (P^) and define the mathematical program {MP): 

k-i ‘^+1 


fc—1 k — 1 

minimize po[y] := (p{xk) + hk'^£{xj,Uj,bj,Xj,Uj,Bj) + E / \\iX„U,,B,)--z{t) 

3=0 3=0 / 


3=0 
bi - bo 


dt 


+ dist^( ,(-oo,M])+dist^(^ ° ,(-oo,M]) 


+ dist 


k-2 

tsll 

3=0 


k-2 


c/,+1 - u, 


(—oo, M]^ + dist^ ^ ^ l^i+i ~ Bj , (—oo, M]^ 


3=0 


subject to equality, inequality, and geometric constraints 


fliv) 

fnv) 

fHv) 

9i{y) 
dp{y) 
y ^ ^ji 
(t>3iy) 


= Xj+i — Xj — hkXj = 0 for j = 0,..., ft — 1, 
= M,+i — Uj — hkUj = 0 for j = 0,..., ft — 1, 
= 6,+i — bj — hkBj = 0 for j = 0,..., ft — 1, 
= {ukt.Xk) - Cfej < 0 for z = 1,..., m, 


X3i\ 


-1 = 0 for j = jr{k),...,f{k), z = l,...,m. 


= {y| 1/2 < ||u,i|| < 3/2} for j = 0,... ,jT-(ft) - 1 and j = j^(ft) + 1,..., ft,, z = l,...,m, 
,, u,j,bj) - z{tj)\\ -e/2<0 for j = 0,..., ft, 


= [Xi,Ui, 


k-l 


3=0 


k-2 


^fc+i(2/):=E / { - z{t) 


</fc+2(2/) := E 


3=0 

k-2 


-U, 


< M + 1, 


4>k+3{y) E |pj + l ~ ^3 < + li 

J=0 

4>k+A{y) ■= Ui-Uo <{M + l)h^, 

(t)k+b{y) ■= bi-bo <{M + l)h'‘, 
y € Ej := {y\ Xj e F{xj,Uj,bj)} for j = 0,... ,fc - 1, 
y G Sfc := {y| xq is fixed, (uo, Co) = (u(0), 5(0))}. 
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It is easy to see that (MP) and (P^) with an arbitrary mapping F are equivalent. Thus y := {z,Z) is 
an optimal solution to (MP), where z := z^ is the solution of (P^) fixed in the theorem, and where Z := 
(Xq, ..., Xk-i, Uq, ..., Bq, ..., Bk-i). Necessary optimality conditions for this type of mathematical 

programs in terms the hrst-order generalized differential constructions constructions of Section 4 are well 
known; see, e.g., [25, Theorem 5.24]. Furthermore, it follows from Theorem 3.4 above that all the inequality 
constraints in (MP) associated with functions 4>j (j = 0,..., fc + 5) are inactive for all k sufficiently large, 
and so the corresponding multipliers do not appear in the optimality conditions. Taking this into account, 
we find A > 0, a e ^ = (^o, • ■ • , a) e Pj = (pj,p“,p5) G as j = 1,..., fc, and 

Vj = (^Oji ■ ■ ■ j ■ ■ ■ I ■ 1 ^{k-l)p^0p ■ ■ • j ^{k-l)j) 

for j = 0,..., fc, which are not all zero and satisfy the conditions (5.6) together with 


(5.10) 


y* G N{y;^j) for j = 0,..., fc. 


m km k—1 

(5.11) -Vo- >^d(po{y) + a,Xg,{y) + EE ^jiVdpiy) + 

i—1 j—0 i—1 j—0 

(5.12) a^gi{y) = 0 for i = l,...,m. 

It follows from the definition of Sj and from fj{y) = 0 that the inclusions in (5.10) are equivalent to 

(5.13) {x*j,uF,bF,X*^) G N(^(xj,Uj,bj, I gP^P) , j = 0,..., fc - 1. 

Note that every other component of j/*, which does not appear in (5.13), is zero. Similarly (a^oj,, Mofc, &ofc) 
can be the only nonzero component of y^. Therefore we have 

(5.14) -y*-...-yl = (-X*o-X*„-xti,...,-4_i,,_i,0,-uSo-«Sfc,---,-4-l.fc-l,0, 

-b*o - b*„ 0, ..., 0,..., 0). 

Let us now calculate the three sums on the right-hand side of (5.11). For notational convenience we just 
specify the nonzero components, which are indexed according to the partition of the vector y introduced 
at the beginning of this proof. This gives us the equalities 




{xk,Uk,bk) 


( k m 

i=0 i=l 

f E 

\i=o 




= (^E > ““j ’ 

2 [Ci? 7 j 0,..., /c, 

-pi if j = 0 

Pj-Pj+i if j = l,...,fc-1 , j = 0, ...,fc, 

Pk if j = k 


k—1 

1,1=0 


= -hkP- 

(X,U,B) 

Introducing the auxiliary Lipschitzian functions (where the first ones are actually smooth around y) 

b+i 


Pjiv) ■= J {Xj,UpBj) - z{t) 


dt, j = 0, ...,fc-l. 
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+ disf 


k-2 

tE| 

3=0 


cr(y) := dist^(^ , (-00, M]^ + dist^ “ 1 ^ ,(-oo,M]^, 


k-2 


c /,+1 - u , 


{— 00 , M]^ + dist^ (E 


\Bj+i - Bj 


3=0 


, {- 00 , M]^ 


and then employing the subdifferential sum rule from [24, Theorem 2.33(c)], we arrive at the inclusion 

k—1 k—1 

dipoiy) C d<f{x'^) + hk'^de{xj,Uj,bj,Xj,Uj,Bj) + ^ Vpj(y) +da{y). 

3=0 j=o 

Since the function dist^(a;; (—00, M]) has the null derivative at all x < M, it implies together with (3.30) 
that da(y) = {0}. Furthermore, the nonzero part of Vpj(y) is given by ^{Xj,Uj,Bj)Pj{y) = '''^ith 

the triple from (5.2), and thus any element of the set Xdipo{y) can be represented in the form 

A( w, d, hkW^, 0, hkw\0, hkV^ + r, + 0“, hkV^ + 9^), 

where f) € dip{x^) and the components of satisfy (5.3). Combining this with the 

gradient expressions above, we deduce from (5.13) via (5.14) the componentwise relationship: 


(5.15) 

-2^00 - 

= XhkW^-pf 


(5.16) 

-^*33 

= XhkW^ + p^ - P^^,, J = 

-1, 

(5.17) 

0 

= Xf) -f y ] CXiUki + Pfc, 

i—1 


(5.18) 

— UqQ — UQk 

= XhkWQ-\-2[^o,uo] -Pi, 


(5.19) 

* 

-^33 

= XhkwJ + 2[^,,m,] -Hp“ -p“+i, j = 

1,... fc 

(5.20) 

0 

= [a,rep^{xk)] +2[^k,Uk] +Pk, 


(5.21) 

^00 

= XhkWQ-p\, 


(5.22) 

-b*- 

^33 

= Xhkw] +p]- p]+i, j = l,...k- 

-1, 

(5.23) 

0 

= -a+pl, 


(5.24) 

~^33 

= X{hkv] +^ 3 )- hkP]+i, j = 0,.. 

. fc — 1, 

(5.25) 

0 

= X{hkvJ + 9J)- /ifep“+i, j = 0,.. 

. fc — 1, 

(5.26) 

0 

= X{hkv] + 9 ))- /ifepj+i, j = 0,... 

,fc- 1. 


Now let us derive from the obtained relationships the necessary optimality conditions of the theorem 
with pq := i^Ok ,UQk,bQk)- We have already got (5.6). Observe now that (5.5) is obviously implied by 
(5.12), the conditions in (5.7) follow from (5.17), (5.20), and (5.23) while those in (5.8) are a consequence 
of (5.25) and (5.26). Arguing by contradiction, suppose that the nontriviality condition (5.4) fails. Then 

it follows from (5.17) that = 0 as well. Since Xq^. = pg = 0, we deduce from (5.15),(5.16), and (5.24) 

that x*j = 0 and X*j = 0 for j = 0,..., /c — 1. Furthermore, (5.25) and (5.26) yield that p" = 0 and 

p^ = 0 for j = 1,..., fc which in turn implies by (5.18), (5.19), (5.21), and (5.22), that also u*j = 0 and 

6*^- = 0 for j = 0,..., fc — 1. As already mentioned, the components of p* different from {x*j , u*j , b*j , X*j ) 
are zero for j = 0,..., fc — 1, and hence y* = 0 for j = 0,..., fc — 1. We similarly conclude that = 0 
due to XQk = Po =0. Getting all this together contradicts the nontriviality conditions in the mathematical 
program [MB) formulated above and thus verifies the claimed nontriviality (5.4). 

It remains to justify of the validity of the discrete-time adjoint conditions in (5.9), which give us a 
discrete-time version of the extended Euler-Lagrange inelusion [23] for the discrete optimal control problems 
under consideration. To get (5.9), we substitute the expressions in (5.15), (5.16), (5.18), (5.19), (5.21), 
(5.22), and (5.24) into the left-hand side of (5.13) and deduce from it that the vector 

{Pj+i -Pj - Ahfcu;J,p“+, - p" - A/ifeu;“,p5+i - pj “ 9^)) - (O, 2 , 0, O) 
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for each j = 0,..., fc — 1 belongs to the normal cone N{{xj, uj , bj, ); gph F). Dividing the obtained 

inclusions by hk > 0, we arrive at (5.9) and thus complete the proof of the theorem. A 

The next theorem directly addresses the discrete approximation problems {P^) for the controlled sweep¬ 
ing process, where the mapping F in (3.20) in given in the particular form (3.1). In this main case of our 
interest we are able to derive, based on the second-order calculations of Section 4, effective necessary 
optimality conditions for (P^) expressed entirely via the problem data. 


Theorem 5.2 (necessary optimality conditions for the discretized sweeping process). In the 

setting of Theorem 5.1 consider the discretized sweeping control problem (Pf), where now the mapping F in 
(3.20) is defined by (3.1). Assume that for the given r.i.l.m. z(-) = (!:(•, u(-), b(-)) of (P'^), which is included 
in (P^), the PLICQ property (3.26) holds and that all the components of u(t) are not zero on [0,T]. Then, 
in addition to the dual elements > 0, = (fo, ■ ■ ■, f,k) G p^ = 

as j = 0,... ,k and subgradients from (5.3) satisfying the relationships in (5.6) 

and (5.8) of Theorem 5.1 with the enhanced nontriviality condition 

(5.27) A'= + ||K'=|| +114'= 11^0, 


there exist vectors pj € R™ as j = 0,..., fc — 1 and € R'” as j = 0,..., fc — 1 such that we have the 
primal and dual/adjoint dynamic relationships 


(5.28) 


^j+1 


m 

i=l 


(5.29) 


xk xk 

Pj + i Pj 

hk 


- X'^wf = 


i=l 


(5.30) 


uk 

Pj+1 


^uk 


hk 


X^wf 



[7^rep„(4)] 




rep„ 


\k(^.xk I ^ axk\ ^xk 

^ ^ ^ ^ Vk^^ ’ ~ 


(5.31) 


bk _ bk 

4 = X^wf - 


hk 


for all j = 0,..., k — 1 together with the implications 


(5.32) 


(5.33) > 0 


< bj, Vji=0 when j = 0,..., fc, i = 1,..., m, 

=» (^vf + - Pj+i^ = 0 j = 0,..., fc - 1, i = 1,..., m. 


(5.34) x’^) < b^j, for all i = 1,... ,m] ^ = 0, j = 0,..., fc - 1, 

as well as the transversality conditions 


(5.35) 


-pt€X>^d^ixt)+Y.pl1ul, 


2=1 


(5.36) 

(5.37) 


pf = - [pl^rep^xD] -2[ek,ut], 


Pki > 0, and {u^i,x^) < ^ = 0 for all i = l,...,m. 
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Proof. In terms of the coderivative construction (4.3) we can rewrite (5.9) in the equivalent form 


(5.38) 


Pj+l 


hk 


-X'^w 


k^,,xk 


^uk 

Pj + l 


hk 


k^,,uk 


X^w 



P% 


Jxk 


hk 


k^,,bk 


X'^w' 






j = 0,... ,fc - 1. 


Using the notation of Theorem 4.3 and employing the coderivative upper estimate (4.13) therein, we deduce 
from (5.38) the existence of vectors r]’- G K™ and 'y^ G QiVj) for which all the relationships in (5.29), (5.30), 
(5.31), (5.33), and (5.34) are satisfied together with the conditions 


— k — k 

(5.39) -(^^)*?y| = —1, where rij € Nm.yy{AjXj -b'^), j = l,...,k-l. 

It obviously follows from the conditions on rjj in (5.39) that G and the implications in (5.32) hold for 
j = 0,..., fc — 1. Defining := ak with ak taken from Theorem 5.1, we deduce from (5.5) that 77 ^ G 
and (5.32) holds for j = k as well. The equations in (5.28) are consequences of those in (5.39) due to the 
definition of Aj given right before Lemma 4.2. The transversality conditions in (5.35)-(5.37) are direct 
consequences of (5.7) due to the relationships in (5.5) and (5.6). 

It remains to verify the enhanced nontriviality condition (5.27). To proceed, suppose that = 0, 
= 0, and pg^ = 0. Then we get from (5.8) that p“*^ = 0 and p^*^ = 0 for all j = 0,..., fc. This ensures 
that yj-' = 0 for j = 0,..., fc — 1 by (5.31), = 0 by (5.35), and consequently pj^' = 0 for j = 0,..., fc 

by (5.29). Furthermore, it follows from (5.30) that =0 for j = 0,..., fc — 1 and from (5.36) that 

Kfei u|] = 0. By recalling the definition of [•, •] in (5.1), the latter conditions readily implies that ^j=0 for 
all J = 0,..., fc. This is due to the assumption of Ui{t) A 0 on [0, T] made in the theorem, which implies 
that A 0 for the discrete approximation due to the uniform convergence of u^{t) —^ u{t) in Theorem 3.4. 
Therefore X^ = 0, = 0, = p^^' = 0 by (5.7), and so pj^ = 0 for j = 0,..., fc — 1. This contradicts the 

nontriviality condition (5.4) in Theorem 5.1 and thus verifies the enhanced one in (5.27). A 


6 Optimality Conditions for the Controlled Sweeping Process 

In this section we proceed with the passage to the limit as fc —00 in the necessary optimality conditions of 
Theorem 5.2 for problems {Pj^) and deriving in this way, with the help of Theorem 3.4 and the developed 
tools of generalized differentiation, necessary optimality conditions in the original optimal control prob- 
lem(s) for the sweeping process formulated in Section 1. Our major case is problem (P^) with 0 < r < T, 
but we also consider the situation when r = 0, which is the same as t = T. The results obtained for {P'^) 
are explicit, i.e., they involve only the problem data and the given local minimizers z(-) while not requiring 
calculations of any auxiliary objects as, e.g., coderivatives. 

Keeping the assumptions above, we impose here some additional ones on the problem data, which 
seem to be reasonable for the controlled sweeping model under consideration and are illustrated below by 
examples. The next theorem uses notation (5.1) together with the symbol ‘co ’ for the convex hull. 

One more remark is needed before the formulation of our main result. Since it is derived by passing 
to the limit in the optimality conditions for the discrete problems (P^), the subdifferential construction 
used in Theorem 5.2 in the case of the nondifferentiable running cost •, •) has to be robust, i.e., outer 
semicontinuous with respect to perturbations of the reference point. As well known, this important property 
holds for our subdifferential (4.4); see, e.g., [24, p. 11]. In the general nonautonomous setting under 
consideration, the robustness of the subdifferential of i with respect to the time parameter is also required, 
and we postulate it in what follows. It does not seem to be restrictive (see the discussion in [23]) and can 
be completely avoided by considering the extended subdijferential of ^ as in [25, Sec. 6.1.5]. 

Theorem 6.1 (nondegenerate necessary optimality conditions for the controlled sweeping pro¬ 
cess). Let z(-) = (a;(-), m(-), 6(-)) be a r.i.l.m. for problem (P^) as t G (0,r) under the assumptions of 
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Theorem 3.4, and let LICQ hold on [0,T]. Suppose in addition that ip is locally Lipschitzian around x(T) 
and the running cost i is represented as 


( 6 . 1 ) 


i{t, z, z) = £i{t, z) + i 2 {±) + u, b), 


where £i, £ 2 , and Ij, are locally Lipschitzian around the given local minimizer z{-) with respect to all but time 
variables. Suppose that ^3 is differentiable in {ii., b) on R” x R with Lipschitz continuous partial derivatives 
and that there is a constant L > 0 such that we have the estimate 


( 6 . 2 ) 


||cJ£ 2 (i, i)|| <for all tG[0,T] and x G R” 

||V4(t,u,6)|| < L(||u|| + || 6 ||) for all t G [0,T] and (m, 6 )gR”+^ 

Suppose finally that z satisfies the estimates 


(6.3) 


lim sup 

k—^oo 


u{hk) - u{0) 


lim sup 

k—^oo 


hk 

u{T) - u(T - hk) 


hk 


<L, 

lim sup 

k—^oo 

<L, 

lim sup 

k—^oc 


bihk)-m 


<L, 


hk 

b{T)-b{T-hk) 


hk 


<L, 


where hk is the step of the uniform discrete mesh taken from (3.4). 

Then there exist a multiplier \>0, an adjoint arc p{-) = {p^, p"^, p^): [0,T] — 
continuous on [0,r] with p{-) G L'^[0,T], signed vector measures 7 = ( 71 ,..., 7 ^) G C*([0,T];R™) and 
^ = (^ 1 ,...,^^) g C*([0,T];R'""), as well as functions {w{-),v{-)) G L“((0, T); R"*) x L2 ((o, T); R") with 


(6.4) 


[w{t),v{t)'^ G cod£(t,z{t),z{t)'^ for a.e. iG[0,T] 


such that the following conditions are satisfied: 

• The PRIMAL-DUAL DYNAMIC RELATIONSHIPS.' 


(6.5) 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 


{ui{t),x{t)) < bi{t) r]i{t) = 0 for a.e. tG[0,r], i = 
rp{t) > 0 ==> (Au“(i) — q^{t), Ui{t)) = 0 for a.e. i G [0, T], i = 1,..., m, 
p{t) = Xw{t) + {0,[-'q{t),\v'^{t) - q^{t)\,0) for a.e. iG[0,T], 

C\ , • 

q^{t) = X—£ 2 (u{t),b{t)), and q’^{t) = X—£ 2 {u{t),b{t)) for a.e. 1G[0,T], 
ou db 


where pf) = (rji {■),... ,r/rn{-)) with the components rii{-) G L^([0, T]; R_|_), i = is a uniquely 

defined vector function determined by the representation 


(6.9) 


x{t) =-'^r]^{t)u^{t) for a.e. iG[0,T], 


and where q: [0,T] jg a function of bounded variation on [0,r] with its left-continuous repre¬ 

sentative given, for all t G [0,r] excepting at most a countable subset, by 

(6.10) q{t) = p{t)-( f ^Ui(s)d 7 *(s), / [rep^(x(s)), ^ 7 ( 5 )] 2 / [u{s),df,{s)],- f dj{s)). 

Jlt.T] Jlt.T] JH.T] ' 


'[t.T] J[t,T\ 

• The TRANSVERSALITY CONDITIONS at the right endpoint: 


( 6 . 11 ) 


-p%T) G A9<^(x(r)) +5]p^(T)i2,(T), 
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( 6 . 12 ) 


i = 1,..., m 


p“(r) +p\{T)x{T) = (^p't{T)+p\{T)x{T)MT))um, 


(6.13) p\T)&W^ and {u^{T),x{T)) <h{T)=^p\{T) = 0, i = 

• The MEASURE NONATOMICITY CONDITIONS.' 

(a) If t G [0, T) and x(t)) < bi(t) for all i = 1,... then there is a neighborhood Vt of t in 

[0, T] such that 'y[V) = 0 for any Borel subset V of Vt- 

(b) Ift G [0,r)U(T —r,r] and 1/2 < ||Mi(t)|| < 3/2 for all i = 1,... ,m, then there is a neighborhood 
Wt of t in [0, r) U (T — t, T] such that f,{W) = 0 for any Borel subset W of Wt- 

• Nontriviality conditions; We always have 

(6.14) A + 11(7(0)11+ ||p(T)||^0. 

Furthermore, the additional assumptions on (ui(0), x(b)} < bi(0) and 1/2 < ||ui(0)|| < 3/2 whenever 
i = 1,... ,m ensure the validity of the enhanced nontriviality condition (X,p{T)) ^ 0. 


Proof. First 'we construct all the functions with the claimed properties satisfying the primal-dual dynamic 
relationships of the theorem. Fix any r G (0,T) and for the given r.i.l.m. z{-) in (P^) consider the discrete 
approximation problems {Pjf) whose optimal solutions = {x^,u^,m) exist by Proposition 3.3 with their 
piecewise linear extensions 0 < t < T, converging to z(-) in the sense of Theorem 3.4. Our aim is 

to derive the claimed necessary optimality conditions for z(-) by passing to the limit from those for ^^(•) 
obtained in Theorem 5.2. To proceed, for each k G N denote by w^{-) and u^(-) the piecewise constant 
extensions to [0,r] of the discrete functions Wj and Vj, respectively, satisfying (5.3). It follows from (5.3) 
and (6.2) due to the strong IF^’^-convergence of the sequence { 2 ^} to z and the local Lipschitz continuity 
of i{t, •, •) that {{w'^(t), v^(t))} is weakly compact in L^([0, T]; =: L^[0,T]. Hence we have 


(w^(-),u^(-)) ^ (w(-),u(-))) weakly in L‘^[0,T] as fc —>• 00 


with some pair (ui(-),u(-)) G L^[0,T]. Employing the aforementioned robustness property of the subdif¬ 
ferential together with the well-known weak convergence result based on Mazur’s theorem (see, e.g., [7, 
Theorem 1.4.1]) allows us to deduce from (5.3) that the convexified inclusion in (6.4) holds. Note also that 
w{-) belongs actually to L°°[0,r] due to its a.e. boundedness on [0,r]. 

Further, based on (5.2) for all fc G N we define the functions 

0xk 

for t G [t], j = 0,..., fc - 1, 

on [0, T] and easily observe by the convexity of the integrand that 


fe-l \laxk\\2 


k-l 


(6.15) 


1=0 1=0 


b+i 


x{t) - 


■'' 1+1 ■''1 




fc-l Ctj+l 


1=0 b 


’i{t) - 


;j.k _ ;j,k 
•^1 + 1 -^1 


hk 


hk 

rT 


dt 


dt = A 


\\x{t) — x^{t)\f dt ^ 0 as fc —1 00 , 


where the convergence is due to Theorem 3.4. This implies that a subsequence of {9^'^{t)} converges to 
zero a.e. on [0, T], The same conclusions hold for the similarly defined functions 0“^(t) and 9^^ ft) on [0, T]. 
It follows from (5.28) that for the piecewise linear interpolations of x^{-) and vf{-) on [0,T] we have 

m 

(6.16) *(t) =for all t G (t^,Ei), j = 0,...,fc-l. 

Now extend lyj’ to [0,T] by (t) := rff for t G [t^,<^_|_i) and define the functions 

(6.17) ’d^ff) := maxjt^l tj <t, 0 < j < k} for all t G [0,T], k gN. 
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We clearly have ^^{t) = for all t G j = 0,... ,k — 1, and 'dfe(f) —t uniformly in [0,T] as 

k oo. The uniform convergence of z^(-) to z(-) on [0,T] readily implies that 

(6.18) ^ z{t) uniformly on [ 0 ,r] as k ^ oo. 

This notation allows us to rewrite (6.16) as 

m 

(6.19) (<) =(t)uf(#(t)) for all t € [0,r]\{tg,...,t^}. 

Consider further the subset of [0, T] given by 

T-.= [0,T]\U 

fcgN 

For any fixed t gT denote by J := I{x{t),u{t),b{t)) the collection of active constraint indices from (3.27) 
and by u^{t) the matrix consisting of the rows it(^('d^(t)), i G J, while u{t) stands for the matrix consisting 
of the rows Ui{t) as i G J. More precisely, this means that uj^(t) = tl 0 (i)(^?^(O) = ^ 4 >(,i){i) for fol 

i € {1,..., I J|}, where | J| signifies the cardinality of J, and where the mapping (j) : {1,..., | J|} —J is a 
bijection. The assumed LICQ condition tells us that the rows of Ui(t) are linearly independent for i G J, 
and consequently we can build the generalized inverse matrix 

u{t) := \u(t)u{t)*'\ ^u{t) for each t gT. 

It follows from (6.18) that uf(d^(t)) — Ui{t) for i G J, and so the generalized inverse 

( 6 . 20 ) u'=(t) := 

is well defined for all k sufficiently large with u^{t) —>• u{t) as fc —>■ oo. We have by the definition of J that 
{ui{t),x{t)) < hi{t) whenever i G {1,... ,m} \ J, and hence (6.18) tells us that 

<b'y{d'^{t)) for f e J" := {l,...,m}\ J 

when k is large. Since d^(t) = tj for some j G {0,..., fc}, we deduce that (^Uj^,Xj) < bj^ for all i G J‘^, 
and so (5.32) yields r]j^ = 0 for this j and all i G J°. Remembering that by construction, we 

conclude that r/f (<) = 0 for all i G J‘^ and large k. This allows us to rewrite (6.19) as 

( 6 . 21 ) = = for large fc, 

iGJ 

where rj'^ (t) collects the components r]f (t) for i G J, i.e., rjf {t) = 77 ^^.^ (t) for i G {1,..., | J|} via the above 
bijection (j). Thus (6.20) and (6.21) ensure the representation 

( 6 . 22 ) {t) = —u^{t)x (t) for large k, 

and the passage to the limit implies that rj^ (t) —>■ —u{t)x (t) as fc —>■ 00 . Define now the required function 
r]{t) = {rii{t),, r]m{t)) for alH S J of full measure on [0, T] by 

(6.23) (t) :=—Ui(t)x (t) for i e {l,..., | J|} and r]i{t):=0 for i G J‘^ 

and observe from the constructions above that 77 ^ (t) —>■ 77 (t) as fc —>■ 00 . Since 77 ^' (t) = 0 for 7 € and 
large k, we get (6.5) by the definition of J and also r]i{t) > 0 for all 7 = 1,..., tti due to r)^ G in 
Theorem 5.2. It follows from (6.21) by passing to the limit that (6.9) holds. The uniqueness of r]{t) in 
(6.9) for a.e. t G [0, T] follows from the imposed LICQ condition. 

To justify the claimed properties of ri{-) in (6.5), it remains to show that ri{-) G L^[0,T]. To see it, let 
us rearrange for each t G T the active components of r](t) by putting rji(t) := 770 (j)(t) for i G {1,..., | J|} 
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and rji{t) = 0 otherwise. Since rji{t) = —Ui{t)x{t) whenever i G | J|} by (6.23) and the convergence 

u^{t) —>■ u{t) is uniform on [0,T] by (6.18), it foiiows that ? 7 ^(-) —>■ ??(•) and hence ? 7 ^'(-) ri{-) strongly in 

i^[0,T]. This not oniy verifies that r]{-) G L^[0,T], but aiso aiiows us to get the estimate 

(6.24) h,J2hj\\^= r\\riHtWdt<M 

i=o do 

with some constant M > 0 independent of k. It immediateiy foiiows from (6.24) and (6.15) that 

(6.25) f ||? 7 f(t)0'^^”(t)|| dt —>• 0 as k ^ oo for ail * = l,...,m. 

Jo 

Next we use the notation of Theorem 5.2 and define q^{-) = {■), {■), {■)) by extending 

piecewise linearly to [0, T] with q'^{t’-) := pj for j = 0,..., k. Construct 7^(-), 'f^(') on [0, T] by 

(6.26) 7 '=(t) := 7 j^, C''(t) := for t G [t^,t^+i), j = 0,..., fc - 1 

with 7 ^(T) := 0 and ^^(T) := Appealing to d^(-) in (6.17), equations (5.29)-(5.31) can be rewritten as 

m 

(6.27) r\t) - 

2=1 


(6.28) qf(t) - X^wf{t) = + 7f (t)x'=(d^'(i)) - vf {t)(t) + - g"'=(d^(t))) 


(6.29) q^^it)-X’‘w^^{t) = --/\t) 

for every t G (t^, t^+i), j = 0,..., fc - 1, and z = 1,..., m, where z?+(t) := for t G , t^+i). 
Now we define p^(-) = [p^^{-)^P^^{-)tP^^{-)^ on [0,T] by setting 

p'=(i):= q\t)+ f (f^^^{s)u’l{§\s)),2[e[s),u\t)\s))] 

(6.30) Jt 

+ [ 7 ns),rep„(x'=(d^s)))], -^\s))ds 
for every t G [0,T]. This gives us p^lT) = q^lT) with the differential relation 

m 

(6.31) /(t) = g"(t) - (5]7.f(0«"(^^t)), 2[e"(t),u"(z?"(t))] + [7^(0, rep„(x'=(d"(t)))], - 7 " 

2=1 

holding for a.e. t G [0,T]. Using this notation, equations (6.27)-(6.29) can be rewritten as 

(6.32) = X'^w^^it), 


(6.33) 


pf{t) = x’^wfit) - v^{t)(^x^v^^{t) + x^e^^{t) - q^'^{d 


(6.34) p^^{t) = X^w’^'^it) 

for every t G (tj,tj+i), j = 0,..., fc — 1, and z = 1,..., m. Define the vector measures by 

(6.35) j d7mes ■= [ 7^(0 [ d^mes ■= I ^^{t)dt for every Borel subset i3c[0,r]. 

Jb Jb Jb Jb 
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From now on we drop for simplicity the index ‘mes’ in the measure notation if no confusion arises. 

Observe next that all the expressions in the statement of Theorem 5.2 are positively homogeneous of 
degree 1 with respect to A^, p^, 7 ^, and Therefore the nontriviality condition (5.27) allows us to 
normalize them by imposing the following relationships whenever fc £ N: 

(6.36) + ||9“'=(0)|| + ||g'’nO)|| + lb'=(r)|| + T 117^^)11^^ + T ll^'=(i)ll* = 1, 

Jo Jo 

which tell us that all the sequential terms in (6.36) are uniformly bounded. Passing below to subsequences 
of fc —?> c» if necessary, we can immediately conclude that A^ — ?► A for some A > 0. Then the equality 
p^{T) = q^{T), the uniform boundedness of the first integral terms in (6.36) and of {w*^)-)} in L°°[0,T] 
implies via (6.27) that the sequence {g^^(-)} has uniformly bounded variations on [0, T], and so it is bounded 
in L°“[0, T], Observe further that the right-hand sides of (6.32) and (6.34) are obviously uniformly bounded 
in L°°[0,T]. Concerning the right-hand sides of (6.33), observe that the sequence {Xp^v^} is uniformly 
bounded in while the remaining summands are uniformly bounded in due to (6.25) and the uniform 
boundedness of Since we also have their uniform integrability by the arguments above, the classical 

Dunford-Pettis theorem on the weak compactness in allows us a subsequence of {_p^(-)}i which weakly 
converges in to some function generating by the Newton-Leibniz formula an absolutely continuous 
function p{-) such that p^{t) —>■ p{t) uniformly on [0,T]. Moreover, the aforementioned Mazur’s theorem 
gives us a subsequence of convex combinations of p^{t) converging to p{t) a.e. pointwise on [ 0 ,T]. 

Using the uniform boundedness of || 7 *’(t)||dt and /J" \\('^(t)\\dt by (6.36), the relationships in (6.27)- 
(6.29) together withp^(T) = q^{T) and (6.25) ensures that {g^(-)} is of uniformly bounded variation. This 
allows us to employ Helly’s selection theorem and find measures 7 £ C'*([0, T]; R™), ^ £ C'*([0, T]; R"*") 
and a function of bounded variation g(-) on [ 0 , T] such that a subsequence of {g^(-)} pointwise converges to 
g(-) while some of {( 7 ^,^^)} weak* converges to ( 7 ,^) in see [3, Definition 3.11, Theorem 3.23, 

and Proposition 3.21]. Thus having q{T) = p{T) and combining it with the a.e. pointwise convergence of 
convex combinations of p^(-) to p{-) justify the possibility of passing to the limit in (6.32)-(6.34) and to 
verify (6.7). Combining (5.33) and (6.15) gives us ( 6 . 6 ) while ( 6 . 8 ) follows from (5.8) and (6.15). 

Now we intend to prove the representation (6.10) for g(-) by passing to the limit in (6.30). It follows 
from the norm convergence of z^{-) z{-) in IU^’^[0,T] that {a;^'(-)} is bounded in L^[0,T] and also that 

supjg[Q_ 7 ^]||'u^'('d^(t)) — it^(t)|| = ^^nmx ^ — M^ll 0 as fc —>■ 00 


with the same for the 6 -components. Furthermore, for any fixed i £ {1,..., m} we have the estimate 

- [ Ui{s)dji{s) 
t Jt 

[ lHs)u'i{^'"{s))ds- [ 7 f(s)u,(s)ds -h [ yf (s)u*(s)ds - [ u,{s)d'y,{s) 

Jt Jt Jt Jt 

f - u,{s))ds + ( Ui{s)^’l{s)ds- [ Ui(s)d'yi(s) , 

Jt Jt Jt 

where the hrst summand vanishes as fc —>■ 00 on [ 0 , T] due to the uniform convergence oi u^{-) to Ui{-) and 
the uniform boundedness of || 7 ^(t)||dt, while the second one vanishes for all t £ [0,T] except for at most 
countably many points because the measures 7 ^' converge weak* to 7 ; see, e.g., [42, p. 325]. Hence 

f 7^(s)uf (i?^(s))(is —> f Ui{s)d'yi{s) as k ^ 00 

Jt Jt 


except for at most countably many t. Proceeding in the same way with the rest of (6.30) gives us (6.10), 
which is the left-continuous representative of the limiting function g(-) of bounded variation on [0 ,T]; see, 
e.g., [3, Theorem 3.8]. Using this and the properties of 'w{-),v{-) and r]{-) established above we complete 
proving all the statements in the primal-dual dynamic relationships of the theorem. 
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Next we justify the transversality conditions, which is a much easier task. Indeed, the validity of (6.11) 
and (6.13) follows by passing to the limit in (5.35) and (5.37), respectively, with taking into account that 
{p^{-)} converges uniformly to p(-) and that p^{T) = q^{T) = Then observe from (5.36) that 

Pk'i + Pfixt = {pI^ + pfixl, ul)ul for all i = 1 ,..., m, 

which gives us ( 6 . 12 ) by passing to the limit as fc —>■ oo. 

Now we proceed with verifying the measure nonatomicity conditions of the theorem. To check the one 
in (a), fix t € [0,T) with {ui{t),x{t)) < bi(t) for i = l,...,m and find a neighborhood Vt of t in [0,T] 
such that for any s G Vt we have {ui{s),x{s)) < bi{s), i = 1,... ,m. This yields {ut {tj), x^ (tj)') < b^{tj), 
i = 1,... ,m, when k is sufficiently large and so € Vj. Thus it follows from (5.34) that for any Borel subset 
V C Vt we have 7 ^(t) = 0 on 17, which implies in turn that ||7mesll(^) = /y ^llTmesll = Iv ll7^(0IM^ = 0- 

Letting k ^ oo shows that ||7||(17) = 0. The measure nonatomicity condition (b) for ^ is justified similarly. 

Our final step is to prove the nontriviality conditions starting with (6.14). Arguing by contradiction, 
suppose that A = 0, g(0) = 0, and p(T) = 0 and hence get A^ —;■ 0, q^{0) -G 0, and p’^{T) = —)• 0 as 

fc —>• oo. Substituting (5.8) into (5.31) with the usage of (5.2), (5.3), and (6.1), we get 


fe-i 


fc-i 


fe-i 


0 A —n j — ^ '' — 

0bk Qbk 
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hk 
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x'^ E IK - Kill + E ^'^IKII + IKII + IKII 
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f^~‘^\\b^ —“Ib^ -Gb^ 

<2AME — K+I+^J- 


\i -0 

k-2 

+ ^E( 

3=0 
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hk 


1 

r - 

1 

r - 

hk 

/ b{t) dt 
Jo 

hk 

/ b{t) dt 

Jr-h’’ 




u^j+2 - 2 m)+i + m }' 

hk 


hk 


A'=E^fcllKll + lKll + lbo 

3=0 


bk\ 


where L is a Lipschitz constant of V£ 2 - Deducing further from (5.8) that 


lbfll<A'=iK^-|l + A^VlKl, 

let us show that ||pi^|| 0 as fc —> oo. Indeed, we get by (5.3) and the smoothness of £2 that = 

^ Combining this with (3.29) and (6.2) ensures the boundedness of {wq^} and hence 

A^IIuq^II ^ 0. It follows from (3.6), (3.29), and (5.2) that is bounded and thus > 0, 

which justifies the claim on ||pi^||. We also have by construction that pg^ = q^(0) -G 0. Involving again 
(3.6) together with (3.30) and the boundedness of by (5.3) gives us || 7 ^(t)||dt —>• 0 as /c —>• 00 . 

Considering now the functions f^(t) from (6.26), we have 


k-1 


ii^nt)iidt = E/^fcf lie 

• /-» ^k 

3=0 


EKi 

3=0 


fc e N. 


It follows from (5.30) due to (3.8) that the estimate 


(6.37) 


Y. iifbi s Y IIp^.i. - II + Y 

3=0 j =0 j =0 

+ E^fei 4 

3=0 


\kaxk 

ykyXk ^ 


Pj+1 


+ E^fcll^Al 

3=0 
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holds whenever i = 1,... ,m. Furthermore, it follows from (5.8) that 


(6.38) 


fc-i 

E IIpS. 

j=o 


■ Pf 




Quk _ Q'^'^ 


hk 


+ X 
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,uk 


uk I 


-<Ml + lbr-Po 


Recalling that ||g^('r)|| = |b^(r)|| —;■ 0 and || 7 ^(t)||dt —>■ 0, we get = 0 on [0,T] by passing to the 
limit in (6.27) due to the weak* convergence of j^(-) in BV]}),T]. Then (6.30) yields p^{-) = 0, and so 

(6.39) max ||(jf^''(t)|| = max ||p“^|| —>• 0 as fc —)• oo. 

te[o,T] j=o,...,k 


Combining (6.39) with (5.2), (6.15), (6.24), and the T^-boundedness of {u“^}, which follows from (6.2) 
and the strong ^-convergence of Theorem 3.4, tells us that the third summand in (6.37) vanishes as 
fc —>• oo. Using the same arguments allowing us to prove that || 7 '"'(t)||(it —>• 0 as fc —> oo, we get by (6.38) 

that the first summand in (6.37) vanishes as well, which therefore verifies that fj' |b^(t)|bt 0. All of 

this leads us to the violation of (6.36) and thus justifies the nontriviality condition (6.14). 

To verify finally the enhanced nontriviality condition under the additional assumptions made, suppose 
by contradiction that (A,p(T)) = 0. By (6.7), p^{t) = 0, p^{t) = 0 for all t G [0,r]. By ( 6 . 8 ), q'^{t) = 0, 
q^{t) = 0 for almost all t G [0,T]. Combining those arguments and (6.10), we get that also q^{t) = 0 for 
almost all t G [0,T]. Using (6.7) again yields that p'^it) = 0 for all t G [0,T]. Therefore, p{t) = 0 for all 
t G [0,T] and q{t) = 0 for almost all t G [0,T]. By using the measure nonatomicity condition, we get also 
< 7 ( 0 ) = 0, hence contradicting the nontriviality condition (6.14) and thns completing the proof. A 

It is worth mentioning (as nsed in Example 7.5 below) that the differentiability assnmption on i with 
respect to (u, 5) can be replaced in the proof of Theorem 6.1 by the following: there is M > 0 snch that 
for all the partitions ^ < ... < tk <T and (u“, u^) G i^£{tj,z{tj),z{tj)), j = 0,... ,k, we have 


(6.40) 


fe-i 


k-l 


Ib “+1 - ^^“11 < M and Y H+I - II ^ 

j=0 j=0 


Indeed, (6.40) is exactly the condition employed above to justify nontriviality (6.14). 


Remark 6.2 (optimality conditions for problem (P)). It is not hard to observe while following 
the limiting procedures developed in Theorems 3.4 and 6.1 that the passage to the limit as r j, 0 in the 
optimality conditions obtained for (P^) in Theorem 6.1 leads us to necessary optimality conditions for 
intermediate local minimizers in problem (P) with the validity of all the relationships (6.4)-(6.14) of this 
theorem hut the second measure nonatomicity condition (b). However, the optimality conditions for (P) 
derived in this way may degenerate in the sense that for any given feasible solution to (P) we can find 
some collection of dnal elements satisfying the nontriviality condition (6.14) such that all the conditions 
(6.4)-(6.14) hold for them. Indeed, this happens when 


A = 0, p{-) = 0, 7 (-) = 0, = (5{o} (Dirac measnre at 0) 

and the adjoint arc q{-) of bounded variation on [ 0 ,T] is constructed as follows: 

_ / (0,-2u(0),0) for t = 0, 
lo fortG(0,r], 


where u(-) is the u-part of the given feasible solution z(-) to (P). Nevertheless, it is important to emphasize 
as illustrated by the examples in Section 7 that, even in the degenerate case, the aforementioned necessary 
optimality conditions allow ns to eliminate nonoptimal solutions and find optimal ones. 


Finally in this section, we consider yet another sweeping optimal control problem much related to (P), 
where the control actions Ui{-) in normal directions are fixed and the optimization is provided by 6 -controls 
changing the position of the moving polyhedron. This problem can be modeled in the following form (P): 

minimize J[x,b] := (p(^x{T)) + J (^£i(t,x{t),b{t),x{t)) + i 2 {b{t))'^dt 
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subject to the constraints in (1.3) and (1.4), where Ui{-) = Ui{-)j i = l,...,m, are fixed absolutely 
continuous functions on [0,r]. Since the equality constraints (1.5) or (1.7) are not imposed, there is 
no difference between problem (P) and its r-perturbations as before. 

We have the following necessary optimality conditions for the new problem under consideration. 


Theorem 6.3 (necessary conditions for problem with fixed normal directions). Let z{-) = 

{x{-),b{-)) be a given r.i.l.m. for problem (P), and let the LICQ condition hold at z{-). Suppose that the 
assumptions of Theorem 6.1 hold whenever appropriate. Then there exist X > 0, an adjoint arc p{-) = 
{p^,p^): [0,r] absolutely continuous on [0,T], L°°-functions (w(-),u(-)) satisfying {w{t),v{t)) G 

cod£{t,z{t),z{t)) for a.e. t G [0,r] with i = £\ +£ 2 , and a measure 7 G C*([0,T];R™) such that for all 
i = 1,... we have the optimality relationships (6.5), ( 6 . 6 ), ( 6 . 8 ), (6.9), (6.11), (6.13), and (6.14) holding 
together with the first measure nonatomicity condition (a). Moreover, (6.7) and (6.10) read as 


p(t) = Xw{t) for a.e. t G [0,T], 

P m p 

q{t)=p{t)-( '^Ui{s)dji{s),- dj{s)) for a.e. t G [0,T]. 

~~i Pt.T] ' 


If finally {ui{0),x{0)) < 6^(0) as i = 1,... ,m, we have the enhanced nontriviality (X,p{T)) 0. 

Proof. Following the proof of Theorem 5.1 shows that (5.2)-(5.8) hold with 0“^ = 0, = 0, = 0 

for all j = 0,..., fc — 1, = 0, and = 0. Inclusion (5.9) reads now as 


'' xk _ xk bk _ bk 

- x^wt. - A'' 


(6.41) 


hk ^ ' hk 

7f.k _ 7j.k . 

-k -k Tk -^i+i 


(7‘ + ^of)) 


G u^, b’j, ^^1^); gphi"(-, ■)), J = 0,..., fc - 1. 


Note that F satisfies the qualification condition in [25, Corollary 3.17], which allows us to deduced that 
D*F,,{x^^,u{t%h)){v) C wo]^^^^mD*F{x),u{t]);b]){v), 

where v := —pj+i + Employing the coderivative estimate (4.13) gives us (5.27)-(5.29), 

(5.31)-(5.35), and (5.37). Then the proof is completed by using the same arguments as in Theorem 6.1. A 


7 Examples and Applications 

We split this section into six examples, which are of a different scale. The first one describes an application 
of the obtained results to a class of elastoplasticity problems, which can be modeled via the sweeping process 
over controlled polyhedral moving sets. The second example addresses a particular sweeping process known 
as the play-and-stop operator, which has various applications to practical models in physics, mechanics, 
engineering, etc. The other examples illustrate special features of the established necessary conditions in 
determining optimal solutions to the controlled sweeping process in one- or two-dimensional settings. 

Example 7.1 (quasistatic elastoplasticity with hardening). We refer the reader to the book [15, 
Chapters 2-4] for models of this type (with no control) and mechanical processes they describe with the 
notation therein; see also some related models in [18]. This example is particularly inspired by models 
in quasistatic small-strain elastoplasticity with hardening. Note that an optimization problem of static 
plasticity with linear kinematic hardening was studied in [17], where the external forces are taken as static 
controls. Here we adopt an essentially different dynamic approach, which seems to be more realistic from 
the viewpoint of mechanical applications. Besides allowing the natural time evolution, we also treat the 
underlying yield criterion as a control action. This leads us to the following optimal control model of 
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dynamic elastoplasticity optimization: to design an elastoplastic material by (dynamically) adjusting its 
yield criterion in order to minimize an appropriate cost. 

To proceed in more detail, consider a body in whose displacement from the initial position is u(t, x). 
The strain e is the symmetric part of the gradient of u, i.e., 

c = i(Vu + V„). 

It can be decomposed into the sum of the plastic strain p and of the elastic strain e by e = e — p. The stress 
(T depends on the elastic strain as cr = Ce, where C is the elasticity tensor and cr satisfies the equilibrium 
equation div cr + d = 0 on an open set with smooth boundary that contains all the possible positions of 
the body together with the boundary condition cr • n = c on dil. Here n denotes the external normal to 
H while d and c represent the external forces that are taken as control actions. We assume that cr = cr(t) 
and p are independent of x. This corresponds to the so-called pseudo-rigid body, see, e.g., [39]. Of course, a 
more realistic model requires dependence on x, but this would lead us to considering the sweeping process 
with an infinite-dimensional state space, which is beyond the scope of this paper. 

If the material undergoes a linear kinematic hardening, then the “plastic flow law” is given by 

(7.1) p € N{(7 — kp; K) 

(see, e.g., [15, pp. 89-90]), where fc is a positive constant, and where K is a compact convex subset of 
called the “region of admissible stresses.” There is a number of interesting practical models of this type 
with a polyhedral region of admissible stresses; e.g., it is a hexagon in the model with the Tresca yield 
criterion described in [15, p. 63]). Denoting q := —kp, we have 

itT := {g G R^l (ct -I- q, Ui) < b^, i = 1,..., m}, 

which induces the polyhedral moving set 

C{t) := {q G R^l {q,u,{t)) < hit) - {a{t),u,{t)), i = 1,.. .,m}. 

This allows us to reformulate model (7.1) as the controlled sweeping process 

-q{t) G N{q{t)-C{t)) 

over the moving polyhedron C{t) with the control functions Ui{t), bi{t), and cr{t). The theory developed 
in this paper can be readily applied to optimize the class of models under consideration with respect to 
general cost functions depending on the state and control variables as well as their velocities. Observe that 
our necessary conditions do not involve d and c directly, but only cr. 

The next example concerns a particular model, which appears in the literature is several contexts. 

Example 7.2 (play-and-stop operator). This name is associated with the sweeping process given by 

(7.2) -xit) G N{xit);b{t) - Z), a:(0) G 6(0) - Z, 

where x G M" (in general x belongs to a Hilbert space), where Z is a closed and convex set (polyhedron in 
our case), and where b: [0,T] —>■ R" is absolutely continuous. We refer the reader to, e.g., [37, Section 7] 
and the bibliographies therein for more details on such operators and their applications. 

To describe the possibility of applying our results, consider for simplicity the case when Z is the 
symmetric rectangle centered at the origin 

Z := {{xi,X 2 ) G R^l \xi\<l3i, \x 2 \ < ^ 2 } with ^i,/32 > 0 

and the control is provided by b{f) = {bi{t), b 2 {t)) for t G [0,T] under the fixed constant u-components 
Ui := (1,0), U 2 '■= (0,1), U 3 := —Ui, U 4 := —M 2 . Then we have 

C{t) = {a; G R^l (a;,Mi) < / 3 i -k 61 (t), {x,U2) < P2 + b2{t), {x,U3) < / 3 i -k6i(t), {x,U 4 ) <132 + b2{t)} 

and are in a position to apply the necessary optimality conditions of Theorem 6.3 to the optimal control 
problem described by (1.2) and (7.2) with the moving set C{t) given above. 
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Now we present several examples illustrating some characteristic features of the necessary optimality 
conditions derived in Theorem 6.1 and also showing how to use these conditions to determine intermediate 
local minimizes. Note that in the examples below the running cost is convex in velocity variables, and so 
there are no difference between intermediate and relaxed intermediate local minimizers. 

The following simple one-dimensional example (as well as the more involved subsequent ones) illustrates 
the procedure of solving problems (P^) whenever 0 < r < T by using Theorem 6.1, even in the case of the 
possible degeneracy for (P) as discussed in Remark 6.2. 

Example 7.3 (calculating optimal controls in one-dimensional problems). Let (P) be given by 

(x — 1)^ • 1 • 

(7.3) n = m = T = 1, xq = 0, <fi{x) := - - -, and £{t, x, u, b, x, ii, b) := -6^. 

It follows from the structure of (P) that we can put u{t) = —1 on [0,1] and thus consider the minimization 
of the cost functional (1.2) with data (7.3) subject to the dynamic constraint 

(7.4) —x{t)€N(^x{t);C{t)^, where C'(t) := {x G R| — x(t) < b(t)} for a.e. tG[0,l]. 

It is easy to see that the variational problem in (7.3) and (7.4) admits an optimal solution; it also follows 
from the general theory due to the convexity and coercivity of the the integrand i in (7.3) with respect to 
velocity. Thus we can apply the necessary optimality conditions of Theorem 6.1 to the problem in (7.3), 

(7.4) and determine in this way its local solution. Employing (6.4)-(6.14) with taking into account that 
(6.12) carries no information in this case give us the following relationships valid for a.e. t G [0,1]: 

(1) w = 0, V = (0,0, b); (2) — x{t) < b{t) => r]{t) = 0; 

(3) r]{t) > 0 => q^{t) = 0; (4) x{t) = 77 ( 1 ); 

(5)(p",p“,p^)(t) = ( 0 ,? 7 ( 09 "(t), 0 ); (6) q-{t) = 0, qHt) = \b{t); 

(7) (g",g“,g'')(t) = (p",p“,/)(<)-(-y' djJ x{s)d^-2j d^,-j dy); 

(8) /(I) >0, -x(l) <&(1)^/(1) = 0; (9)-p"(l) = A(x(l)-l) - /(I); 

(10)A+||g(0)|| + |b(l)||,^0. 

We consider first that case where x{t) ^ 0 for a.e. t G [0,1]. Then it is evident that x{t) = —6(1) > 0 
and hence it follows from (4) that q^(t') = 0 for a.e. t. Observe furthermore that (5) implies that p(-) is 
constant, which ensures by (7) that dy is constant on (0,T] as well. The latter means that either y 
is zero or it is a Dirac measure concentrated at 1 = 1. In both cases we have that g^(l) is constant by (5) 
and (7), and so b{-) is constant by ( 6 ) provided that A 7 ^ 0; otherwise we do not have enough information 
to proceed. Assuming A 0 yields in this case that there is only one feasible trajectory satisfying the 
necessary optimality conditions; namely {x{t),u{t),b{t)) = (1/2, —1, —1/2) with the cost value of 1/4. The 
case of A = 0, where no information can be deduced on x, cannot be ruled out. To examine finally the 
opposite case of x(-) = b{-) = 0, we see by the same arguments as above that q{-) is constant on (0, T] with 
q^{t) = A — p^(l). This choice satisfies necessary optimality conditions with the cost value of 1/2 > 1/4, 
and thus we found a reasonable candidate to be an optimal solution to this problem. 

Observe that the problem in Example 7.3 can be also treated by necessary optimality conditions of 
the Pontryagin Maximum Principle from conventional control theory by taking into account that the state 
constraint therein is active, i.e., x(l) = —6(1), and so we can consider 6 ( 1 ) as the new control. In fact, PMP 
allows us to show in this setting that (x(l), ?2(1), 6(1)) = (1/2, —1, —1/2) is a global minimizer. However, it 
is not the case in the following modification of the previous example, where the moving constraint is not 
active, and we cannot reduce the sweeping process to a conventional control system. 

Example 7.4 (necessary conditions for the controlled sweeping process versus PMP). The only 
difference of this example in comparison with Example 7.3 is that the running cost is given now by 

£(l,x,u, 6 ,x,u, 6 ) := ^(( 6 - 1 )^ -f | 6 p). 
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The trivial choice of b{t) = 1 for all t S [0,r], gives us the value 1/2 of the cost function (1.2). If instead 
the moving constraint is active, then obviously u{t) = —1 and x{t) = —h{t) > 0 on [ 0 , T], which shows that 


(^( 1 ) - 1 )^ 
2 




It is easy to see that the trivial solution above satisfies the necessary conditions in Theorem 6.1 (take 
g^(-) = p^{-) = A = 1 and let the other dual elements vanish). Since in this case the sweeping state 
constraint is not active, we cannot employ PMP as in Example 7.3. 


The next example shows how to exclude nonoptimal solutions by using necessary optimality conditions 
from Theorem 6.3. Observe that in this example the measure 7 has an atom and the corresponding adjoint 
arc q{-) is discontinuous inside the time interval. 


Example 7.5 (excluding nonoptimal solutions). Consider problem (P) with u = —1 on [0,1], n, to, T 
as in (7.3), xq = 1/5, <p{x) := (x — 1)^, and the running cost given by 


£(t,x,b,x,b):=(b + t — So(t))^ + ajb + 4t — 2j with o > 0, so(t) := < 


(‘py 


i( 


t- lift + 


iff < i, 

if i i 

otherwise. 


It is easy to see that the couple {x{t),b{t)) = 
a = 0 , where 


vo{t) := 


[vo{t),—t + so{t)) on [0,T] is optimal for this problem if 

if^<i 

It ifi<t<l, 

I I otherwise. 


Let us show it is not (at least in the intermediate local sense) if a > 0. It can done it by applying the 
necessary optimality conditions of Theorem 6.1 while observing that all its assumptions are satisfied with 
the validity of (6.40) for £2 nondifferentiable in b. Furthermore, due to a;(0)it = —1/5 < 1/25 = 6(0) we 
can employ the enhanced nontriviality condition (A,p(l)) ^ 0. 

To proceed, deduce from (6.9) and ( 6 . 6 ) that q^{t) = 0 for a.e. t £ [1/5,4/5]. Moreover, since r] = 0 
both in [0,1/5] and in [4/5,1] thanks to (6.5), we deduce from (6.7) that p{-) = 0, and so p{t) = p{\) for 
all t € [0,1]. Assuming by contradiction that A ^ 0, we get from ( 6 . 8 ) that q^{t) = aA for t > 3/4 and 
q^{t) = —aX for t < 3/4. Then (6.10) tells us that 


q^{t)= I d 7 (s) and q^{t) = — f d'y{s) for a.e. te[ 0 ,l]. 

Therefore, on one hand the measure 7 is zero on (1/5,4/5), while on the other it must have a nonzero 
mass at t = 3/4. This is a contradiction, which shows that A = 0. at the same time we get p(l) = 0 by 
the transversality conditions (6.11)-(6.13) due to ^(l)!! = —4/5 < —1 + 6/25 = 6(1). This contradicts the 
enhanced nontriviality and thus verifies that (x(-), 6 (-)) is not optimal for (P) if a > 0. 


Finally, we present a two-dimensional example that can be analyzed on the basis of Theorem 6.3. 
Example 7.6 (controlled sweeping process in two dimensions). Let the data of (P) be: 

lIxlP • 1 • • 

(7.5) n = TO = 2 , a:o = ( 1 ,1), T = 1, (/5(x) = —^, and £(t, x, m, 6 , ±, ft, 6 ) :=-( 6 ^-b 62 ) • 

Consider the version (P) of this problem with the fixed normal vectors ui = (1,0) and U 2 = (0,1) and 
apply Theorem 6.3 to determine optimal solutions b{t) = (61 (t), 62 ( 0 ) x{t) = (xi(t), X 2 (t)) on [0,1]. 
The necessary optimality conditions of Theorem 6.3 give us the relationships on [0,1]: 
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(1) w{-) = 0, v-{-) = 0, = {hi-),ki-)); (2) ^ 0 gf(i) = 0, z = 1,2. 

(3) p^(-) is constant with nonnegative components, and —pf (•) = Xx{l)+p^{-)ui, z = 1, 2 is also constant. 

(4) g"(t)=p"- 7 ([t,l]), r 7 '’(<) = A 6 (t)=p' + 7 ([t,l])fora.e. te[0,T]. 

(5) A+ || 9 ( 0 )|| + ||p(l)|| ^ 0 with A > 0. 


Observe first that the trivial solution with x{t) = (1,1) and b{t) = 0 on [0,1] satisfies necessary conditions 
(take Pi = P 2 = ~lj p\ = p\ = h = l 2 = 0, and A = 1). In this case the cost value is 1. If the i-th 
constraint is pushing (i.e., Xi{t) < 0 on a set of positive measure), it follows from (4) that 7 i([f, 1]) is 
constant on that set and also bi is constant for i = 1,2 provided that A 7 ^ 0, which is supposed to hold. 
We consider only the simplified case where pushing occurs on at most one interval. There are the following 
three possibilities in this case: 

(a) Both constraints are pushing with constant speed at the same time in the interval [0,i?], where 

0 < < 1 is to be determined. 

(b) The constraints are pushing alternatively (with constant speed); by symmetry we may assume that 

they push for the same time, say first ui in the interval [ 0 ,zl] and then U 2 in the interval [i?, 2 i?], where 
0 < < 1/2 is to be determined. 

(c) Only one constraint is pushing in the interval [0,i?], again with constant speed; by symmetry we 
may assume that the first one is active. 

To proceed further, denote the constant speed of the f-th moving constraint by /3i < 0. Then in case 
(a) the cost value is calculated by 


i (^2 + ^)(/ 32 +; 3 |) + ^(/3^+/32) + l. 


and it is subject to minimization over / 3 i ,/32 < 0 with 0 < zl < 1. Straightforward calculations show that 


/3i = 132 


-1 

1 + 1 ?’ 


and in this case the i?-component of the gradient of the cost function is negative, and so i? = 1. Thus 
/3i = P 2 = which gives us the cost value A. 

In case (b) we have the same cost value as in (a) while 1 ? G (0,1/2]. The same calculations tell us that 
the optimal cost with this strategy is obtained for 1 ? = 5 and /3i = /32 = — 5 and its value is In case (c) 
the cost value is | with the choice of i? = 1, /3i = — ^, and /32 = 0. Combining all the above allows us to 
conclude that the strategy in case (a) is the most appropriate when A 7 A 0. Finally, observe that if A = 0, 
we do not have enough information to proceed. 
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